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Abstract 

In this paper, we describe Bohr-Sommerfeld rules for semi-classical 
completely integrable systems with 2 degrees of freedom with non degen- 
erate singularities (Morse-Bott singularities) under the assumption that 
the energy level of the first Hamiltonian is non singular. The more singu- 
lar case of focus-focus singularities is studied in |27| and (26). The case 
of 1 degree of freedom has been studied in |lo| . 

Our theory is applied to some famous examples: the geodesies of the 
ellipsoid, the 1 : 2-resonance, and Schrodinger operators on the sphere 
S 2 . A numerical test shows that the semi-classical Bohr-Sommerfeld rules 
match very accurately the "purely quantum" computations. 
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Introduction 

In this paper, we describe extensions of results of jlCj to completely integrable 
semi-classical systems with 2 degrees of freedom. If H± , Hi are two commuting 
/i-pseudo-differential operators on a 2-manifold A, we introduce the momentum 
map F — (Hi,H 2 ) : T*X — > R 2 where Hj (j = 1, 2) is the principal symbol of 
Hj] we assume that F is a proper map. If o = (0,0) is not a critical value of 
H, the existence and the construction of solutions of the system HjU = 0(h°°) 
has a long history and existence of solutions is well known to depend on the 
Bohr-Sommerfeld rules involving action integrals and Maslov indices of loops 
generating the homology of the fibre of F which is a 2-torus. 

We will assume that o is a critical value and more precisely that the critical 
points are of Morse-Bott type. A very simple classification of such points with 
the corresponding normal forms is given in ]25[ ] . The focus-focus case have been 
described in p^j. We will be interested in the case where the set of critical 
points of H is a 2-dimensional manifold with a transversally hyperbolic (saddle) 
singularity: it means that is not a critical value of H\ and that H2 restricted to 
Poincare sections of the flow of Hi admits critical points of saddle type. The set 
of critical values of F is then a 1-dimensional submanifold of R 2 . The main result 
of our paper is a description of singular Bohr-Sommerfeld rules in this situation. 
These rules give necessary and sufficient conditions for existence of solutions 
of the system HjU = 0(h°°), j = 1,2 and approximations of the solutions. 
More precisely, we show the existence near the singular fibre A D = _F _1 (o) of a 
Hamiltonian H p with periodic flow which allows to reduce the classical and also 
the semi-classical study to the 1-dimensional case on the reduced phase space. 
This reduced phase space can be singular because the S 1 action induced by the 
flow of Hp is not principal in general; non trivial isotropy group isomorphic to 
Z/2Z may appear. 

We provide a description of A G : the topological type can be rather unusual 
like a Klein bottle! 

The precise description of the quantisation rules is given in theorem 2.7. 
There is one rule giving a quantum number associated to the periodic orbits 
of Hp and another rule given in terms of the graph G which is the quotient of 
A Q by the S^-action associated with H p . With a viewpoint similar to |27| but 
somewhat different from p0| , we interpret these rules as a universal regularisa- 
tion of the usual Bohr-Sommerfeld rules for tori. This, in addition, allows us 
to have a description of the joint spectrum inside a fixed neighbourhood of o, 
but raises some technical difficulties due to the possible non-connectedness of 
the fibres of F, especially in the C°° category. 

At the end we describe 3 examples for which we provide explicit calculation 
and numerical checking: 

1. High energy limit for eigenvalues of the Laplacian on ellipsoids. 

2. Semi-excited states for anharmonic oscillators with a resonance 2:1. 

3. High energy limit for the Schrodinger spectrum on the 2-sphere. 

For the last two examples, numerical computations of eigenvalues of large ma- 
trices are compared with the eigenvalues obtained from the singular Bohr- 
Sommerfeld rules. We observe a very good accuracy of the results even for 
not very big quantum numbers. 
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1 Classics 

The goal of this section is to give a semi-global description of the Lagrangian 
fibres of a 2-degrees-of-freedom integrable system having non-degenerate rank- 
one singularities of hyperbolic type. 

Let (M, oj) be a symplectic manifold of dimension 4, and let Hi, H 2 be two 
Poisson commuting Hamiltonian functions in C°° (M, R) . The corresponding 
momentum map will be denoted by F = (Hi, H2); we shall always assume F to 
be proper, Hi to be non-singular on the level set A D := F~ 1 (o) (for some point 
o € R 2 ). Moreover, we assume that A G is connected and that the critical points 
of F on A are transversally non-degenerate, in the following sense: 

for any Poincare section E of Xi , the restriction of H2 to E is a Morse 
function. 

The main results of this section are 



The description of the topology of the fibre A Q (section 1.2) 



The construction of partial action-angle coordinates in a full neighbour- 
hood of A Q ; we show in particular that there exists a Hamiltonian H p de- 
fined in some neighbourhood of A that Poisson commutes with Hj , j — 



1,2, and all orbits of which are periodic (theorem 1.6). Moreover, up to 
finite covering, a neighbourhood of A is symplectomorphic to a product 
of T*S 1 by a "global" Poincare section. 

The construction of normal forms for the system near each connected 



component of the critical set of F (theorem 1.13) 



The core of these results can be obtained from Fomenko's classification |Q, 
and the first two points are actually due to Nguyen Tien Dung p3| . How- 
ever, we felt that these results deserved a detailed description with independent 
arguments. 



1.1 Notations 

For convenience of the reader we group here some of the notations: our in- 
tegrable system is given by the momentum map F = (H\,H2) ■ M — > R 2 . 
Because F is proper, it is a momentum map for a Hamiltonian action of R 2 on 
M. o — (a ,b ) is a critical value of F, U small disk around o in R 2 . The 
smooth energy level is 5 Q = H^ l (a a ). Q, is F -1 (J7); A Q = F _1 (o). 7 = LlJI^j 
is the critical set of F in A Q . T = UTi is the critical set in fL A Q \ 7 = UA^ .-, 
where ^L^x is a smooth Lagrangian cylinder whose closure is A^ - x U 7; U jj. 
For any Hamiltonian H somet hing we will denote by X SO mething the associated 
Hamiltonian vector field. G is a graph with TV vertices associated to A Q . The 
Hessian of H2 restricted to E is 7is(i?2)- The absolute value of its determinant 
(with respect to the density induced by the symplectic for on E) is independent 
on E and denoted by |7Ys (^2) | • 



1.2 Topology of A 

Proposition 1.1 The critical set 7 
dimension 1 of A Q which is a finite 



of F in S a is a compact submanifold of 
union of disjoint periodic orbits 7, (i — 
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1, ■ • • , N) of X\. The ji admits orbit-cylinders I\ which consists of Ji(a), a 
close to a a . 

Proof. Locally the reduced manifold of S := {Hi = a } is symplcctomorphic 
to a Poincare section E of X\. Because X\ is transversal to E in S , a neigh- 
bourhood V of E in S is diffeomorphic to E x /, for a small interval /, in such 
a way that the trajectories of X\ are of the form {a} x I, a e E. Because H 2 
is constant under the A^-flow, /^({c} x I) = (H2)fs(cr). By hypothesis, H 2 
restricted to E admits a non degenerated critical point at a :— 7 n E. Since 
the rank of F if invariant by the flow of Xi, we must have 7 n V = {a} x I, 
which says that 7 is a smooth manifold of dimension 1. Because, by definition, 
7 is closed in the compact set A D and hence is compact, it must be a union of 
circles 7 = [Jji. Then each of these circles must be an Ai-orbit, and only a 
finite number N of such critical circles arises due to the properness of F. 

The non-degeneracy of (i^fs implies that the isolated critical point a = 
<r(a ) depends smoothly on a close to a Q . Therefore the above description 
extends to any leaf of the foliation Hi — a, a close to a , yielding a smooth 
family of circles 7i(a). Since a small neighbourhood of S in M is diffeomorphic 
to S x M such that a is a coordinate for the second factor, the union I\; = 
U a close to a a 7i( a ) is diffeomorphic to a cylinder S 1 x (R, a ). □ 




Hi 



Figure f : Image of the momentum map. The set of critical values near o is a 
smooth curve parameterised by Hi. If the critical point is a saddle point, the 
regular values are on both sides of the curve of critical values. If it is a maximum 
or a minimum, only one side is occupied. 

If the critical point a is a local maximum or minimum of H2, A Q reduces 
to one elliptic periodic orbit, a situation which was studied a long time ago by 
several people (see 0). We will from now assume that this critical point is a 
saddle point. Then one can have several critical circles in 7; we show now how 
they can be connected to each other inside A„. 

Proposition 1.2 A D \7 is a union ofM. 2 -orbits that are cylinders. Each cylinder 
contains in its closure 1 or 2 ji's. We will denote by = Ufc{^{ij}} ^ e 

set of cylinders that connect 7, and 7j . 

Proof. Since 7 is M 2 -invariant, A D \ 7 is a union of orbits, on all of which the 
action is non-singular. These orbits are therefore 2-dimensional quotients of R 2 , 
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hence tori, cylinders or planes. Because A Q is connected and contains singular 
points, any orbit contains critical points in its closure, which excludes tori. The 
local structure near each 7, will show independently the existence of periodic 
sub-orbits, leading to the non-triviality of the stabilisers of the R 2 -action, which 
excludes planes. Hence all orbits are cylinders. 

Now, the closure of such a cylinder consists of singular orbits: there are 1 or 
2 of them. □ 

Definition 1.1 We define the graph G as follows: G has N vertices (where N 
is the number of critical circles in A a ), and there are exactly \Aii t j\\ different 
edges connecting the vertices i and j . 

Let £ be a Poincare section at a point m e 7^. By hypothesis, £ n A D is 
diffeomorphic to a "hyperbolic cross", the union the local stable and unstable 
manifolds W ± (m) for the flow of Xjj 2 ■ Let SI be a small neighbourhood of 7,, 

and define W (ji) as the union of the connected components of (A \t;) fl fi 
intersecting W ± (m). These manifolds do not depend on the choice of m. 

Proposition 1.3 Either {^i) andW~(~/i) are diffeomorphic to the disjoint 
union of 2 cylinders or both are diffeomorphic to 1 cylinder. In the first case 
the vertex i of G has degree 4 while in the second it has degree 2. 




Figure 2: The vertices of G are of degree 2 or 4. 

Proof. The 2-manifold W + (ji) = W + (ji) U 7.; is a bundle on 7; whose fibre is 
an interval. There are exactly 2 possibilities up to diffeomorphism: the trivial 
and the Moebius bundle. In the first case, removing 7^ gives 2 cylinders, while in 
the second it gives only 1 cylinder. Both bundle are isomorphic because the sum 
of their tangent bundle on 7^ is a R 2 — bundle that is orientable as a symplectic 
bundle. □ 

Definition 1.2 In the first case, 7* is called direct, in the second case, it is 
called reverse. 
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a) direct case b) reverse case 

Figure 3: The neighbourhood of a critical circle 



1.3 The classical commutant 

Lemma 1.4 For any smooth function K commuting with Hi and H2, the vector 
field Xk can be uniquely written (in a neighbourhood of Ao ) 

%k = a A"i + b A2 , 



for smooth functions a, b commuting with Hi and i?2- 



Proof . Near any nonsingular point of F, apply the Darboux-Caratheodory 
theorem which makes Hi = £ and H2 = 77 in a local symplectic coordinate chart 
(2:, y, £, 77). Then K = K(^,n), and the result follows. 

Near a critical point in some 7^, we use Theorem L5 below that reduces the 
situation to Hi = £ + a Q and H2 — yn), and hence to (Hi = £, H2 = yrf), 
for which we use the result of 25, Lemme 2.6] (p(| Lemme 2.2.7]). □ 



Theorem 1.5 Around any point m £ 7,, there exists a canonical chart {x, y, £, rf) 
in which 

H\-a a = £ 

for some smooth function <I> defined near the origin, with 

<9 2 $(0,0)>0. (1) 



Proof. First construct a Darboux-Caratheodory chart for Hi, i.e. a system 
of symplectic coordinates (x, y, £, rf), with canonical form loq, in which £ = 
Hi — a a . In these coordinates, the plane {0} x {0} x K 2 is a Poincare section 
for Hi, hence by hypothesis (y,n) — > H2(S,,y,i]) has, for each small £, a non- 
degenerate saddle point. The application of the Isochore Morse Lemma [jll| 
with parameter £ yields a local symplectomorphism 4>^ (y, rf) of the (y, ?7)-space, 
depending smoothly on £, such that 

H 2 (x,^^(y,r,)) = ^,yr,), (2) 

for a function $ with 82^(0, 0) ^ 0. Applying the canonical transformation 
(Vi v) ~^ (~V> v) if necessary, one can assume that $ satisfies the condition (|l|) of 
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the theorem. The map (j> '■ (%, £> U, v) l— * ( x ) 0? (2A 7 7)) is a local diffeomorphism 
but need not be symplectic. A modification of x shall solve the problem. 
The 2-form w\ := 4>*(u>o) splits as follows : 

uii = uio + d£ A /?, 

for some 1-form /3 = (3(£,y,rf). From dwj = 0, one gets d£ A d(3 = 0, which 
means dr ViV \@ = 0. Hence there exists a smooth function f(£,y,T)) such that 
d(y,ri)f = ^" ^ e ^ ne diffeomorphism : 



so that 



<f>i{x,t,v,y) ■= - f(^y,v),^y,v), (3) 



^wi = 0Jwo + 0*i (d£ A 0) = #wo + d£ A /? (4) 



= LU -d£Adf + d£Af3 = LU - (5) 

Thus 0O0! is symplectic, and because it does not change i?i = £ — a D , it answers 
the question. □ 



1.4 A periodic Hamiltonian H p in f2 

Theorem 1.6 ([p3|) There exists a unique (up to additive constant) Hamilto- 
nian Hp in n that fulfils the following conditions: 

1. Hp Poisson commutes with Hi and Hi. 

2. The flow of Hp is 2ir periodic with minimal period lix outside T. 

3. OnT t , 

X p = a i (H 1 )X 1 , (6) 

with oti a positive function. 

Then the flow of H p is 2ir-periodic onTi if the vertex {i} is direct, and n -periodic 
in the reverse case. 

Lemma 1.7 Near each ji, there exists a unique (up to some additive constant 
if a Q =/=0) Hamiltonian Hii = Hi — \i(H\)H\ which is critical on IV 

Proof. If dHi = fi(Hi)dHi on r<, one gets the differential equation 

t\' t + X t = fi(t) 

which always admits a local solution. □ 

Definition 1.3 A class of path z — » [y z ] (where [ ] means a homotopy or homol- 
ogy class) is called smooth if there exists locally a representative that smoothly 
depends on z. 
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Lemma 1.8 Let us denote by L z the orbit of z by the R 2 — action. There exists 
a unique mapping z — > [y z ] that is a smooth map from 51 into Hi(L Zl 1) and 
such that, if z £ ji, [y z ] — f [7i] where ji is oriented according X\ and v = 1 
(resp. 2) if is direct (resp. reverse). 

Proof, a) The main point for this proof is to construct s uch a smooth family 



[7 Z ] near 7^ By the Morse-Bott lemma (see Appendix 1.7) applied to the 



Hamiltonian iJ 2 ,i of Lemma 1.7, there are a coordinate X\ on the circle 7$ = 
(R/Z) and a fibre bundle F± of dimension 2 on 7, defined as quotient of the 
trivial bundle on R by identifying (xi,w) € ji x R 2 with (x\ + 1, ±w) such that 
-£?2,i = 2:3^4. We can assume that X\ has the same orientation as d/dxi on 7.;. 

The above Morse-Bott lemma can be applied with the parameter Hi . There- 
fore, one gets coordinates (x±, X2, £3, X4) on a full neighbourhood of ji in M by 
letting x-2 = H\. 

We choose then j z to be the path given in these coordinates by t — > (t, X2, £3, £4) 
with t € R/Z (direct case) and t E R/2Z (reverse case). These path are drawn 
on the leaf (Hi,H2.i) = const, and hence on a leaf of F = (iii,^)- The last 



assumption of lemma 1.8 follows from proposition 1.3 



b) Far from 7,, we construct 7^ on A Q by taking the cylinder equator with 
the orientation which in the affine structure of is given by projecting 
X 1 . Then we extend to the nearby Lagrangian leaves by local triviality of the 
foliation by orbits. Because these j z are homotopic to the ones constructed in 
a) and lying on the same Lagrangian leaf, it is easy to realise this homotopy as 
an isotopy, thus yielding a smooth family of loops in f2. □ 

c) It remains now to use these j z in order to define the action variable. 

Lemma 1.9 The symplectic form uj is exact in 17, i.e. there exists a, 1-form 
in £1 with da = uj. 

Proof. The set A Q is Lagrangian and any 2-cycle can be deformed inside A Q . 
□ 

Proof (of Theorem Put H p (z) = (1/2tt) / a. H p is smooth and 

commutes with Hi and H% (because it is constant on R 2 — orbits). Moreover on 
R 2 — orbits that are tori the orbits of X p are 2ir— periodic with orbits homotopic 
to 7 Z (by usual action-angle coordinates). On the 7,'s, the period is 2ir in the 
direct case and ir in the reverse case. 

The affine structure on the Lagrangian cylinders ^ and the condition (3.) 
implies the uniqueness of X p on A Q . Now suppose H' is another Hamiltonian 
with the same properties, and let z be a point in A Q \ 7. Since the orbits under 
X p and Xp of z are equal, the orbits of points near z in a same level set of F 
(different from A D ) are homotopic. But these level sets are Liouville tori for 
which we know that X p and X' p must be equal. □ 

Remark 1.1. Step a) of the proof does not use the nature of Ar, 1 . Therefore, 



a)+c) gives a Hamiltonian H p verifying the conditions of Theorem L6 but in 



a neighbourhood of 7^ only. This suffices to prove that the R 2 -action has non- 
trivial stabi liser s, whence A/ 4 .-. must be cylinders, thus finishing the proof of 



proposition 1.2. A 
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1.5 S 1 reduction 

The flow of X p yields a locally free Hamiltonian action of S 1 on f2, which is free 
outside r. Let c a = H p (m) = H p (A ). We denote by W the reduced space 

W = i/ J 7 1 (c Q ) n Q/S 1 . 

W is a symplectic orbifold (see eg. It is a smooth manifold if and only 

if the action is free, that is if and only if no vertex of reverse type are present 
in the graph G. Otherwise, it has singularities at the critical orbits 7^. Since 
these critical orbits come in families depending on the value of H±, yielding 
local orbit cylinders, and because dH p (m) = XdHi(m), A 7^ 0, only one orbit of 
each local cylinder meets ff~ 1 (c ), ensuring that the critical orbits give isolated 
singularities in W. 
Let 

H q = -b(H 1 -a ) + a(H 2 -b ), (7) 

where a > and b are the real constants such that, on A G , X p = aX\ + bXi 
(cf. Lemma |L1|). Then A D = H~ x (c ) n H' 1 ^). (This still holds of course 
for a generic choice of (a, b).) Since H q is ^-invariant, it defines a smooth 
Hamiltonian function H q on W. The graph G can be viewed as the quotient of 
A a by S , and thus is identified to the level set i?~ 1 (0). 

Proposition 1.10 If the S 1 -action is free (i.e. all vertices of G of degree 
4 ), then a neighbourhood Q of A a in M is diffeomorphic to the direct prod- 
uct S 1 x R x W (hence A Q is diffeomorphic to the direct product S 1 x G - 
these diffeomorphisms are equivariant with respect to the natural action of S 1 
on itself). 

Remark 1.2. In this case, W can be regarded as a "global" Poincare section 
for X p . A 

Proof . We choose now ft to be of the form o x / where I is some open 
interval around 0, fl is a small invariant neighbourhood of A G in H p 1 (c ), and 
H p (Vt x {£}) — c a = £. If the action is free, S7 is a principal S^-bundle over W. 
It is topologically classified by its holonomy class in H 1 (W, Si), where S% is the 
sheaf of germs of smooth functions on W with values in S 1 (see (l9|). Using the 
short exact sequence O^Z^R^S^^O, and the fact that the sheaf Ro is 
fine, one gets an isomorphism 

yielding the so-called Chern class of the bundle. But W retracts onto G and G 
is 1-dimensional, so H 2 {W,1) = 0, and tt is a trivial bundle. □ 



Theorem 1.11 ( [ 23 1 ) If the S 1 - action is free, then (il,u>) is symplectomorph 



ic 



to a neighbourhood of S 1 x {0} x W in 

{T*S 1 xW, d4Adx + 7r*uj w ), 

with H p — c Q = £. Here 7r is the projection onto W and low is the symplectic 
form of W . 
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Proof. First apply proposition 1.1C and let £ = H p — c a be a coordinate for the 
K factor. Then choose the conjugate angle variable x (pick up some coordinate 
9 in S , an origin 9q, and let x(9) be the time required to go from 9q to 9 under 
the Hamiltonian action of £), so that {£,2;} = 1. Because u> is ^-invariant, it 
does not depend on x\ using the equivariant Darboux-Weinstein theorem p8| |, 
one can assume that u> — u)\^—q, and so does not depend on £ either. 

Because for any £, X v is w-orthogonal to S 1 x {£} x W, one easily checks by 
taking local coordinates on W that 

lo = d£ A dx + d£ A ir* f3 + tt*u>w, 

where is a one-form on W, and u>w a non-degenerate 2-form on W. The 
closedness of u (and its independence on £) implies dwyy = and df3 = 0, the 
latter yielding d^ir* (3) = d£ A 7r*/3. Let us now apply Moser's path method to 
get rid of this term. We let 

u>t := g?£ A dx + ir*u>w + td(^ir* /?) , 

and wish to construct an isotopy Lp t of diffcomorphisms of Q such that ifttot = 
ujq. (ft is then given as the flow of the vector field X t defined by ix t 0Jt+£,^* P — 0. 
It is easy to check that u>t is non-degenerate for all t so that X t is uniquely 
defined. Moreover, because of its defining equation, X t is of the form £t*Yt (l 
is the inclusion W ^ T*S 1 x W), where Y t is a vector field on W satisfying 
iY t (u>t) \W + P — 0. Therefore ip t is of the form 

{x,£,w) (x,^(j>^ t (w)), 

(where <f>t is the flow of Y t ) and hence preserves x and £. If Yt can be integrated 
up to the time to > 0, then X t can be integrated up to the time 1 for £ ^ to. 
The diffeomorphism <£>i then answers the question. □ 

Remark 1.3. The formula u) = d£Adx + ir*uJw ensures that ojw is the natural 
symplectic form on W obtained by the reduction process. A 

Theorem 1.12 In the general case, there exists a smooth double covering fl* 
of f2 in which the action is free. The reduced manifold W* is a covering of W 
that is ramified of degree 2 at the critical orbits ji . 

Proof. Choose VI to be a relatively compact invariant neighbourhood of A D in 
M. Then H~ 1 {c o )C\0, has a smooth non-empty invariant boundary, and because 

this boundary does not meet any critical orbit, the closure W = 7? ? T 1 (c ) D VI/ S 1 
is a relatively compact surface with a non-empty smooth boundary. Let pi, i = 
1 ... I be the images under reduction of the critical circles 7* and let W be the 
surface W after removal of small disks Di around each pi. It is still a smooth 
surface with boundary, whose fundamental group is free Q, and generated by 
some n\,...,Hk,8\,...,8i, with Si = dDi. Let V C n\{W) be the free subgroup 
generated by m, . . . , fik, 8\, . . . , <5f , and let W* be the corresponding smooth 
covering of W . Gluing back the disks Di defines a covering space W* of W that 
is ramified of degree two at each pi. Define now fT C fl x W* such that the 
following diagram commutes: 

O* W* 



Q W 
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Since the local structure near 7, (see Theorem 1.13 in the next section) gives 
a model for the covering Q* — > Q, Q* is naturally endowed with a smooth 
structure compatible with that of Q. The lifts of 7^ are critical circles in 51* 
that become of direct type. As a result, all critical circles in 57* are of direct 
type, which means that the lifted S 1 action is free. □ 



1.6 Normal forms near 7, 

Choose any Hamiltonian function H q near ji that commutes with Hi and H2 
and such that, for any o in the local curve of critical values of F, 

• H q (A ) = 0; 

• on A Q \ 7j, X p and X q are linearly independent, and the automorphism: 
[X\^X-i) — > (X p ,X q ) is orientation preserving. 

For instance, the previously defined H q (eq. ((7J)) is a good choice, which is 
independent on i, but any generic linear combination H q = a{H\ — a ) + ji(Hi — 
b ), (3 > would also do. 

Then the following theorem states a simple normal form near 7, for the new 
system (H p ,H q ). It will be the main tool for the semi-classical analysis, for it 
reduces the situation to the case of a cotangent bundle. 

Theorem 1.13 There exists coordinates (x, y) on R = Rx]A, B[ that give local 
coordinates on W + (~fi) by taking the quotient by 

(x,y) -> (x + 27r,y) 

in the direct case and 

(x,y) -> (x + ir,-y) 

in the reverse case and a canonical diffeomorphism of a neighbourhood of '7^ into 
a neighbourhood of the % = c a -section" ofT*(W + (ji)) (recall that c a = H p (A ) ) 
such that with respect to canonical coordinates we get 

/ H p = £ 

\ H q = myr,), 

for some smooth function $ defined near the origin, with <92 < I > (0,0) > 0. 

Proof. 1. We first wish to prove that the restriction of H q to the locally reduced 
manifold W := H p ~ 1 (c ) / exp tX p has a non-degenerate saddle point. Let m € 7; 
and let (s, a, u, v) be local coordinates near m such that a = H\ — a and the flow 
of X\ is just translation on the s variable. Then because dH p (m) = \dHi{ni) 
with A ^ 0, the map (s, cr, u, v) — > (s, H p — c Q , u, v) is a local diffeomorphism of 
(M 4 , 0) that sends (s, 0, 0, 0) to (s, 0, 0, 0). Therefore one can take (u, v) as local 
coordinates for W, and we wish to prove that (H q ) \w(u, v) has a non-degenerate 
saddle point at (0,0). 



Let M = 




be the matrix of smooth functions such that 



{dH u dH 2 ) = M ■ (dH p , dH q ). 
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Since H2 is of the form H2 = F(Hi,u,v), one has 

dH 2 = K.dHx + A, 

where K = d\F and A = A(H\,u, v) is a one-form on {0} x {0} x R 2 depending 
smoothly on H±, that vanishes at m and whose differential at m is a non- 
degenerate quadratic form of hyperbolic type. Using M. one gets 

(c - Ka)dH p = (Kb - d)dH q + A. 

The claim is that (Kb — d)(m) ^ 0. Indeed, (c — Ka) and (Kb — d) cannot 
simultaneously vanish because M. is invertible. It suffices then to see that 
A(m) — whereas by hypothesis dH p (m) ^ 0. 
Therefore, we have, on TW C kerdH p , 

W lTW = -^- d A, 

which implies that d(H q )^w possesses, as A does, a non-degenerate differential 
of hyperbolic type. 

2. Now, we consider a neighbourhood of the whole critical circle 7$ and use 
Weinstein's theorem with S 1 action to reduce to the cotangent space of W + (^i) 
with (Hp — c = £,,H q = y, rj)). In these coordinates (we consider first 

the direct case), W is naturally identified with {0} x {0} x IR 2 , so the previous 
point shows that (y,rj) — > H q (^,y,r]) has, for each small £, a non-degenerate 



saddle point. Then the proof goes exactly as that of Theorem 1.5, Equations 
(||) through (H) are valid if Hi and H2 are replaced by H p and H q , and loq is 
the canonical 2-form of T*(W + (-fi)). 

In the reverse case, the proof is the same but we need the Isochore Morse 
lemma for functions that are invariant by the involution <j(x) = —x: it it then 
possible to choose the diffcomorphism F commuting with a. This fact follows 
easily from the proof given in jflj: in the lemme principal p. 283, we choose rj 
such that <J*(n) = —rj. It implies that X t commutes with the involution. □ 

Remark 1.4. We decide to give to the graph G the orientation of the flow of 
H q defined by (0). Near a vertex 7^ , it is also given by the flow of the normal 
form yr\. A 



1.7 Appendix: Morse-Bott lemma 

Definition 1.4 Let f : X — > R be a smooth function. A submanifold W of X 
is called a Morse-Bott critical manifold if every point w S W is a critical point 
of f and if the restriction of f ' (w) to the normal bundle T W X /T W W is non 
degenerate. 

Morse-Bott critical manifold arises in many situations especially when / is 
invariant by a Lie group action. An extension of the Morse lemma is available 
in that case. In some situations, there is global topological problem with the 
subbundles N± of the normal bundle generated by eigenspaces of /" associated 
with > (resp. < 0) eigenvalues. 
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Lemma 1.14 (Morse-Bott lemma) Assume we have a Morse-Bott connected 
critical manifold W for a function f : X — > WL. Let N be the normal bundle 
of W and F the Hessian of f which is a non degenerate quadratic form on N. 
Then there exists a diffeomorphism of a neighbourhood ofW on a neighbourhood 
of the 0-section in N which conjugates f to F + c. 

If W is connected, a complete set of invariant of /, up to smooth conjugacy 
near W, is given by the pair (N + , N_) of bundles on W up to isomorphism. 
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2 Semi-classics 

The aim of this section is to express the singular Bohr-Sommerfeld quantisa- 
tion rules for quantum integrable systems whose classical counterpart fulfils the 
hypothesis of the previous section. 

Let X be a 2-dimensional differential manifold, and let H 1 (h) 1 H 2 (h) be 
commuting /i-pseudo-differential operators, with real principal symbols Hi, H 2 . 
Assume that the momentum map F — (Hi, H 2 ) satisfie s the hypothesis of sec- 
tion [j]. In all of this section, the 1-form a of Lemma L9 is taken to be the 



can onical Liouvillc 1-form of the cotangent bundle T*X. Then H p (Theorem 
|l.6|) is uniquely defined as the action integral with respect to a. For any E in 
the image of F, the sub-principal form ke is the closed diffcrentiable 1-form on 
Ae ■= F^ 1 (E) defined at its regular points by ke(Xj) = ~ r ji where rj is the 
sub-principal symbol of Hj. 

2.1 The microlocal normal form 



We will prove here a semi-classical analogue of Theorem 1.13] , which was partic- 



ularly fit for this purpose since it reduced the situation to that of a cotangent 
bundle, for which the usual pseudo-differential quantisation can be used. A 



semi-classical analogue of Theorem 1.11 should also be interesting, but would 



involve symplectically reduced cotangent bundles, for which Toeplitz quantisa- 
tion is needed, a theory that we don't want to enter here. 



In this section a critical circle jj C A Q is fixed. Theorem 1.13 identifies, via 
a symplectomorphism tp, a neighbourhood of 7^ in T*X with a neighbourhood 
of the zero section of a cotangent bundle of the form T*(M 2 /cr), where a(x,y) — 
(x + 2ir, y) in the direct case, and <j(x, y) = (x + it, —y) in the reverse case. 

It is easy to check that Weyl quantisation satisfies, for a symbol a € Cfi°(M. 2 ) : 

a*0 P f(a) = OpY(aoT*cx)a*, 

where a* is the adjoint operator u ^ uo a, and T*a is the cotangent lift of a. 
Therefore, if a = a o T*er, then OpF (a) acts on the space of functions u that 
are invariant under a : u o a = u, which is the space of functions defined on 
a cylinder in the direct case, and on the Mcebius strip in the reverse case. In 
particular, Qi(h) — Op^ (£) and Q 2 (h) — Op]^ (yrj) are well-defined differential 
operators on M. 2 /a : 

w*)-75- «*> = T("i + 2)- M 

Let ^° be the algebra of operators of the form Op^ r (a(h)) for classical sym- 
bols a(h) on T*(R 2 /er), modulo those whose symbol is 0(h°°). For more details 
about the microlocal calculus and 0(h°°) remainders, see eg. ||, [p6| . Before 
stating the result of this section, we introduce the following spaces : 

Definition 2.1 The classical and semi-classical commutants <L c i(lj) an d£h("/j) 
are defined as follows : 

={/ G C°°(r*(M 2 /a)), {/,£} = {f, m } = near^}; 
£h(lj) = {P(h) S *°, [P, Qi] and [P, Q 2 ] are 0(h°°) near 7,}; 
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(Recall that in our coordinates, jj — ((R x {0})/er) x {c } x {0}, where c = 
-ffp(A ).) Because the symbols of Qi and Q2 are polynomials of degree 2, the 
operators in &h{lj) are exactly the Weyl quantisations of symbols of the form 
J2 h k a k with ak € <£ c /(7j')- 

Theorem 2.1 (Microlocal normal form) There exists an elliptic Fourier in- 
tegral operator U(h) associated to the canonical transformation ip of Theorem 



1.1 S , an invertible 2x2 matrix M{h) of pseudo- differential operators in <Lh{lj), 
and complex-valued functions of h : t\{h) and e 2 (h) admitting an asymptotic 
expansion in C[/i]: 



e 1 (h) ~ x: <P ^)~E4 fl 

smc/i i/iat, microlocally near 7j : 

ET 1 ^ - a OJ ff 2 - 6 )t7 = M.(Qi - e x ,Q 2 - he 2 ) + 0{h°°). (9) 

If Hi and H 2 are formally self-adjoint, then U (h) can be chosen to be microlo- 
cally unitary, and the functions tj are real-valued. 

The first terms of e\(K) (of order respectively h and h 1 ) are given by the 
formulae 

^=c = ^J s a; (10) 

e<p = ±- [ « o + /i(<0/4, (11) 

where fi is the Maslov index of any regular part of A a , and S is any cycle 
associated to an S 1 -orbit on A a \ 7 (and recall that k is the sub-principal 
form of the system). 



The first term of e 2 {h) is given by the formula 

(0) _ ( M - r 2 
\Hv{H 2 )\V* 



(12) 

where A is defined in Lemma \l.1{ ( recall that ri is the sub-principal symbol 
of Hi andH.s:(H 2 ) is the transversal Hessian ofH 2 ). Note that Tlj:(H 2 ) is 
also equal to the (y,rj)- Hessian of H 2 j (the latter was defined along with 
A in Lemma 

Remark 2.1. Recall that there is a choice of sign in the canonical chart ip of 



Theorem 1.13 . If the other sign is chosen, then e 2 becomes —e 2 . A 

Proof. Consider the direct case first. The proof goes along the same lines as 
in HH, so we skip the details and indicate only the crucial points. First take U 
as any Fourier integral operator associated to ip (note that by construction this 
symplectomorphism is exact in the sense that it preserv e the action integral). 
Since Hi and H 2 commute with H p and H q , Theorem 1.13 implies that the 



principal symbols of U 1 H\U and U 1 H 2 U are in the classical commutant 
^ci{lj)- The following division lemma is easily proved as in p5| : 
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Lemma 2.2 Any function K £ € c i(jj) that vanish on jj can be written (in a 
neighbourhood of"fj) : 

K(x, y, rf) = K(£, y, rf) = a{£ - c c ) + byn, 

for some smooth functions a and b in £ c /(7j). 

Applying this lemma to Hi o jp and H 2 ° tp solves the principal part of equation 

The next steps are obtained by conjugating U by elliptic pseudo-differential 
operators, yielding transport equations of the form : 

Lemma 2.3 Given any functions (r\,r2) such that 

{ri,yrj} = { r 2,£,}, 
there exists K\ , K 2 £ £ c i (7j ) arid a function f such that : 

/} = K l - r i and {yv, f} = K 2 -r 2 . 

Proof. Let 

i r 2 * 

Ki {£,■> Vi v) = 7r~ I n ( x , f > y, v)dx. 
27r Jo 

Of course, {£, K\} = 0, and using the hypothesis of the lemma, one has 



1 f 2n 

{Ki,yrj} = — l {r 2 ,£}dx 
27T Jo 



2vr J dx 

Now, let 

fi(x,€,y,v) = / K i(ty,v) - n(x',C,y,v)dx'. 

Jo 

Then f\ is a smooth function on T*(R 2 /er) that satisfies - using the same kind 
of calculation as above : 

{fi,yv}(x,£,,y,v) = r 2 (x,£,,y,n) - r 2 (0, f, y, rf). 

Then the wanted function / is sought under the form f — fi + fi , which leads 
to the system 



{£,/ 2 } = 0, and 



{Wl,h} = K 2 - r 2 (0,£,2/,T?). 

It suffices to see £ as a parameter and apply a known lemma in the (y, rf)- 
variables (see eg. Theorem 2]). □ 
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In the reverse case, the proof of the theorem is the same provided we deal with 
functions that are invariant under T*o~. But if K (£, —y, —if) = K(£, y, 77), lemma 



2.2 still applies, yielding functions a and b with the same properties. The same 
is true for the transport equation. Then each step of the proof can be quantised 
via Weyl's formula to yield well-defined pseudo-differential operators on M 2 /a. 
Thus the result still holds for the reverse case. 

The proofs for formulas ([h]), (JTl|) and ([l2]) are given in section 2.4 , but the 



formula for ej is apparent in the proof of Theorem |2.4| below, and the formula 
for 62 can be directly checked using the fact that the subprincipal symbol is 
preserved under conjugation by an elliptic Fourier integral operator at a critical 
point of the principal symbol. □ 

2.2 Microlocal solutions 

We investigate here the solutions of the system 

(Hi{h)-a )u = 0(h°°), (H 2 {h)-b )u = 0(h°°), (13) 

microlocally on a neighbourhood of the critical Lagrangian A Q . If the operators 
Hj depend smoothly on some additional parameter E € R d that leaves the 
principal symbols intact, then all the results presented here depend smoothly 
and locally uniformly on E. This applies in particular to the investigation of 
the joint spectrum in a window of size O(h) around (a Q , b ), where Hj is to be 
replaced by H — hEj . 



Theorem 2.1, applied to all critical circles 7^, yields a finite set of semi- 
classical invariants (ei, 7 (/i), £2, , (/&))■ We show here how these quantities are 
related to the solutions of fll3|). 

Theorem 2.4 (The global quantum number) 

• The asymptotic series e±(h) = eijQi), modulo KL, depend neither on j, 
nor on the particular way to achieve the normal form of Theorem \2.A. 



• The system ( \19[ ) admits a microlocal solution near any (and then all) 
S 1 -orbit (including critical circles) if and only if the following condition 
holds : 

ei (h) €hZ + 0(h°°). (14) 
Remark 2.2. Since e\(h) is determined by Hi(h) and H2(h), the fulfilment of 



equation (14) seems to impose a quantisation condition on h. While we can stick 
here to this interpretation, another possibility would be to recall that everything 
(and in particular e\(h)) smoothly depends on the point o in the curve of critical 
values of F. Then equation ( p"4j ) can be interpreted as a quantisation condition 
on o, which leaves h free to vary in a full neighbourhood of 0. This viewpoint 



is made clear in section |2.4| (cf. Corollary p. 14 ) . A 



Proof of Theorem 2.4. Following [|26|, we introduce the sheaf (£, A Q ) of germs 
of microlocal solutions on A D , as a sheaf of C^-modules, where is the ring of 
all complex functions of h, c(h), such that 



\c(h)\ < C.h 



-N 
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for some constants C,N, modulo those functions that are 0(h°°). Note that 
the vector operator F = {Hi^Hz) acts on the huge sheaf over A G of germs of 
all admissible distributions modulo microlocal equivalence, and (£, A D ) can be 
seen as the kernel of F. The question is to find out how local germs can be 
glue d together to form a nontrivial global section of (£, A G ), i.e. a solution of 
near A . 

It was shown in |^6| that the restriction (£, A G \ 7) to the non-singular points 
of F is a locally constant sheaf, and the germs £(p) at any non-singular point 
p form a free module of rank 1, generated by a standard WKB solution. The 
existence of nontrivial global sections of (£, A Q \ 7) is then characterised by 
the nullity (mod. 0(h°°)) of the associated holonomy (or "Bohr-Sommerfeld 
cocycle") : 

X(h) € H 1 (A \7,R/2ttZ). 

Since A G \ 7 is a disjoint union of cylinders A| whose homology fli(Af i3 -j) 
is generated by the cycle represented by any oriented S' 1 -orbit, we get a finite 
set of holonomies ^(h) characterising (£, A D \7). 

Apply now Theorem |2.l[ The system (Q) is then, on a neighbourhood of 
7j, equivalent to the following standard system : 

Qiu = e± t jU, Q 2 u = he 2 ju. (15) 

At any non-singular point p£ Q \ jj , the standard WKB u solution generating 
£(p) is therefore of the form 

u(j>)=e i ° J t Lx v(y), (x,y) ~p£ R 2 /<t. 

This implies that 

— [ X(h) = ^4^- + 0{h°°) (modZ), 
2tt J s h 

where 5 is the cycle on A D \ 7 associated with the orbit S 1 (p) . This proves 

1. that ei,j(h) does not depend on the particular way to achieve the normal 
form; 

2. that J s X(h) does not depend on the choice of the Lagrangian cylinder 
containing p - i.e. for j fixed, all X% j\ (h) are equal - which in turn 
proves 

3. that ei — e\ t j does not depend on the choice of the critical circle 7^ (since 
A Q is connected). 

The first part of the theorem is now proved. Moreover, the condition ([TiJ) is 
necessary and sufficient for the existence of a (non-trivial) solution near a regular 
orbit. This condition remains therefore necessary for the existence of a solution 
near a critical circle. We are thus left with the proof of the sufficiency of this 
condition for critical circles, which is achieved by the next proposition. □ 

Proposition 2.5 Let 7^ C 7 be a critical circle. Let d ~ 2 or 4 be its degree in 
G(A D ). If the condition (jl^) is fulfilled, then the set £(7?) of germs of microlocal 
solutions on jj is a free Ch-module of rank ^. 
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Proof. Let n = n(h) £ Z be such that t\ — hn + 0(h°°), and let p £ 7 ? -. We 
know from ^ proposition 17] that the module of microlocal solutions of ( |l5| ) at 
p is free of rank 2, generated by 

u± d = e m * [ l ±B>0 _Le^ ln W j . (16) 

V vlyl / 



If 7j is direct, this immediately implies that the module £(7j) °f microlocal 
solutions of ([l3]) on the whole circle 7, is also free and of rank 2. 

In the reverse case, the distribution + C_tt_ on T x R is invariant 

under the involution it if and only if it has the parity of n in the variable y, 
which reads here 

c _ = e in^ C+ _ 

£(7j) is in this case a free module of rank 1, and its generator depends on the 
parity of n. □ 



2.3 The abstract Bohr-Sommerfeld rules 

We assume here that the first condition ( [l4| ) is fulfilled, and show that the 
existence of global solutions of ( p"5| ) can be read on the graph G = G(A Q ). As 
before, let n — n(h) £ Z be such that t\ = n + 0(h°°), 



Because of Theorem 2.4, for any point p £ A Q , there exists a microlocal 
solution on a neighbourhood of the orbit S 1 (p). We shall use this fact to con- 
struct from the sheaf (£, A Q ) a new sheaf (£, G) on G C W (recall that W is 



the symplectic orbifold of section 1.5) that will in some way encode whether 
(£, A e ) has a global section. To each point p £ G we associate the free module 
£(p) generated by the germs of standard basis at p, which will be of rank 1, 
as follows. For details of the construction when all vertices are of direct type, 
see |0). 

Denote by the vertex of G corresponding to the orbit jj , and let 7 = 1J 7^ . 

• At a regular orbit in A D \ 7, a standard basis is just any basis of the space 
of solutions near jj, so we let £(p) = £,(S 1 (p)). 

• At a vertex 7, of degree 4, we use the definition of a standard basis from 
[ [To! , which we recall here. 

The edges connecting 7^ are oriented according to the flow of yrj. More- 
over, near jj, W is a smooth oriented surface (the orientation i s giv en by 
minus the symplectic form). It is shown in H that Proposition 2.5 in the 



direct case can be restated as follows : let I1I2 (resp. -Z3/4) be the disjoint 
union of the two local edges leaving 7^ (resp. arriving at 7^) with cyclic 
order (1,3,2,4) - with respect to the orientation of W near 7^. Z(Iil2) 
and £(13/4) are free modules of rank 2. Then there exists a linear map 
Tj : £(13/4) — » £(11/2) such that u is a solution in a neighbourhood of jj 
if and only if its restrictions satisfy u^r^j^ — TjU^n^y In other words, 
if we "feed" the system with two functions on the entering edges ^3 and 
I4, then these functions are propagated on the leaving edges I\ and I2 in 
a unique way (see fig.0). 
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4 3 2 1 

Figure 4: Propagation of solutions at vertices of degree 4 and 2 



Tj ; as a 2 x 2 matrix 



, 4, and express 
(defined modulo 0(/i°°)). Moreover, one 



One can choose a basis element for each £(-£,-), i = 1, 
a 6 

c G? 

can show that the entries are all non-vanishing. It is then easy to check 
that a new choice for the basis elements does not change the cross-ratio 
pj = j£. In our situation, pj can be explicitly calculated; one finds (§) : 



p j = p j (h) = -e 2 ^ h \ 



(17) 



The choice of a matrix Tj fixes the basis elements up to their multiplication 
by a same factor. We shall call the choice of the basis elements of £(/»), 
i = 1, . . . , 4, a standard basis whenever Tj has the following expression : 



'2nh 



r(/3)e 



Plnh 



-i0% 
J0§ 



with /3 = h + if2,j- or with the notations of fief] 



with 



T 3 =T(e 2J (h)) 



£j = S(e 2 ,j(h)) 



1 



3 i/3f 



1 



(18) 



(19) 



(20) 



Remark 2.3. In (T(J the factor e -1 ^ was absent in the definition of Tj 
( p^[ ) . Its introduction here will greatly simplify the treatment of Maslov 
indices (see also ||). A 

Remark 2.4. Equation ( |l7|) proves that e 2 j is a semi-classical invariant 
(modulo iZ) of the critical circle jj : it does not depend on the particular 
way the normal form is achieved. A 



Remark 2.5. As it is presented here, the notion of a standard basis 
seems to be attached to the graph G endowed with a specific labelling at 
vertices of degree 4. The form of the matrix Tj shows that the different 
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possible labellings of the four hyperbolic branches i = 1,2,3,4 yield 
the same set of standard basis, provided I\ and I2 are the local unstable 
manifolds (for the flow of yrj - which means for the flow of H q ), I3 and 
Ii are the local stable manifolds, and on the oriented manifold W (which 
is smooth at vertices of direct type) the branches appear in cyclic order 
(1,3,2,4). Furthermore, it can be easily checked using the standard ba- 
sis (S) given below and the fact that T\f = f (where /(y) = f(—y)) 

that T{e)T{— e) = ( ? ! ]. Therefore, exchanging the local un/stable 



x 1 / 

manifolds just amounts to changing the sign of e% j. A 

• At a vertex 7,- of degree 2, the space of solutions has dimension 1, so 
we could just, as in the regular case, call any solution a standard basis. 
However, in order to isolate the "singular" components of the holonomy, 
we prefer the following convention which is more in accordance with the 
previous case (degree 4). 

Let I u and I s be the local unstable and stable manifolds of jj . A choice of 
basis elements (e u , e s ) for £(I U ) and £(/ s ) will be called a standard basis 
if 

with 

Cj = C{n,e 2 ^(h))^e-^e- m ^£ 3 (l + i(-l) n e- e2 - 7T ). (22) 
-|r( / 3)e /?ln,l cos(|( / 3 + n)). 
Notice that Cj depends on n mod 4. 
Remark 2.6. If £2 j € M., one has 

g. _ \ e iargr(l/2+ic2.j)+ie2,i hi(h) (23) 

3 Vi + 

Therefore Tj is unitary and \Cj\ = 1. A 

(£, G) is a locally flat sheaf of rank-one modules, and hence is characterised 
by its holonomy: 

hol:Hi(G)->C h . 

In terms of Cech cohomology, if 7 is a loop in G, and fii, O2, ... , fig — Qi is an 
ordered sequence of open sets covering the image of 7, each Oj being equipped 
with a standard basis Ui, then 

hol(7) d = £1,22:2,3 ■ ■ ■ Xt-i,l, 



where dehned in by m = %i.jUj on Qi n £1 



3 ' 



Definition 2.2 (The singular invariants) 1. The "principal value" k a of 
the sub-principal form k is the cocycle on A Q defined as follows ( see also 

my- 
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if [A, B] C A D \ 7 is a non singular path, then 



[A,B] 



[A,B] 



• if [A, B] C A Q is a path intersecting once and transversally a unique 
critical circle jj, and oriented according to the flow of H q (ie. A is 
on the local stable manifold and B is on the local unstable manifold) 
then 



[A,B] 



lim 

a, b — >m 



[A.a] 



[b,B] 



•4°> 



where R a ,b is the parallelogram (defined in any coordinate system) 
built on the vectors ma and mS (m — [A, B] D 7,- - see fig.^l. 




2. The "regularized" Maslov index fi on A„ is defined as follows: 

• The contribution from a regular path in A a \ 7 is the usual Maslov 
index of the path. 

• Let S = [A, B\ C A D is a small path intersecting once and transver- 
sally a unique critical circle 7/, such that A belongs to one of the 
hyperbolic branches (1,3,2,4) and B to an adjacent branch (ie. 5 
makes a turn of angle ±§J- S can be continuously deformed into a 
path St drawn on a regular leaf of F. Then the usual Maslov index 
for this path is constant for t small enough (So = 5), and we define 

A*(«5)-.= M(*t)±Q + X{*=2}nJ, (24) 

where ± = "+ " if 5 turns in the direct sense ( with respect to the cyclic 
order of the branches) and "—" otherwise, and X{ C ; J . = 2} = 1 */ 7j is 
of degree 2 and X{dj=2} — otherwise. 



2 SEMI-CLASSICS 



23 




Figure 6: The local Maslov correction (for a vertex of degree 2, replace f by 



Remark 2.7. The sign in (gj) is negative if 5 is oriented according to the 
flow of H q and 5t belongs to a region in phase space where H q > 0, and changes 
whenever one of these two conditions changes. Of course the rule of Fig. ^ is 
simpler to use, but this correspondence will be used for the proofs of section 



2.4. A 



Proposition 2.6 1. The holonomy of the sheaf (£, G) has the form hoi = 
e l l s ( h )}/ h ^ where [n*9(h)] € H 1 (A D , (C, +)) admits an asymptotic expansion 
in non-negative powers of h. (it : A D — > G is the projection associated to 
the S 1 -reduction.) 

2. Let ^2g> [6i]h e be the asymptotic expansion of [0(h)] := n*[6(h)]. Then 
the first two terms are given by the following formulae : 



[#o] = [ct] (the Liouville 1-form on A ); 

• = [£<,] + M|- 

Proof. We just prove here the existence of the claimed asymptotic expansion. 



The formulae for 9 will follow from our refined analysis in Section 2.4 (Corollary 



2.21) 



For this purpose, it suffices to show that one can choose local sections u a of 
(£, G) for which the transition constants c a ,p have the required form. On the 
edges of G, this follows from the regular theory of WKB solutions (see eg.]27|]). 



At a critical circle, we apply the normal form (Theorem 2.1), and choose the 
following standard basis (u± is defined in Eq.(|l6|)): 



in the direct case, 



/2 



F -r {u -^ (25) 



t2\ 



• in the reverse case, 



e u = e% + (-l) n e 2 
e s = i n e 3 + (-i)"e 4 . 
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We see then that the restrictions of these solutions to any edge are standard 
WKB solutions, whose phases admit an asymptotic expansion in powers of h. 
□ 

The dimension bi = dim IT (G) is given by Euler-Poincare formula : 

bi = #{edges of G} - N + 1 

(recall that N is the number of vertices of G) . Moreover, if we write N = p + q 
with p the number of tetravalent vertices and q the number of divalent vertices, 
then it is easy to see that 

#{edges of G} = 2p + q, 

so that &i = p+1. 

We can now cut b\ edges of G, each one corresponding to a cycle Si in a 
basis {5i, . . . ,5^) of H 1 (G), in such a way that the remaining graph is a tree T 
(Hi(T) = 0). Then the sheaf (£, T) has a nontrivial global section, i.e. there 
exists a standard basis u a on each edge e a such that they extend to a standard 
basis at each vertex. 

Theorem 2.7 (£, A Q ) has a nontrivial global section if and only if the following 
linear system of 3p + q + 1 equations with the Zp + q + 1 unknowns (x a € 
< C-h)ae{edges of T} has a nontrivial solution : 

1. if the edges (a\, a%, a>2, ct±) connect at a tetravalent vertex jj (with the 
prescribed orientation), then 

{x a3 , x a4 ) — Tj {x ai , x a2 ) , 

2. if the edges (a u ,a s ) connect at a divalent vertex (with the prescribed 
orientation), then 

Xa B — CjX au ; 

3. if a a and ai are the extremities of a cut cycle Si, then 

x ao = hol(Si)x ai . 

Here we assume the following orientation: Si can be represented by a closed 
path starting on the edge a n and ending on a\. 

Remark 2.8. When solving the system, it is immediate (if p ^ 0) to replace 
the equations of type 2. and 3. into those of type 1., in order to finally obtain 
a linear system of size (2p) x (2p). If p = then equations of type 2. and 3. 
combine together to yield a unique equation in one variable. A 

Proof . Any global section u of (£, A D ) can be characterised by the set of 
constants x a € C/j such that Wf£( eQ ) = x a u a . By definition of the standard 
basis (u a ), conditions (1.) and (2.) are necessary and sufficient for (x a u a ) to 
extend to a solution near every critical circle jj. It remains to check that the 
solutions at the extremities ao, a\ of a cut cycle Si can be consistently glued back 
together. Since (u a ) is a global section of (£, T), u ai is the parallel transport 
of u ao along Si, which means that, as local sections of (£, G) (or (£, A Q )), they 
satisfy u ai = ho\(Si)u ao . Therefore the solutions x ao u ao and x ai u ai can be 
glued back if and only if condition (3.) holds. □ 
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2.4 The spectral problem 

The goal of this section is to investigate uniform estimates for our system when 
it depends on spectral parameters. Specifically, we look now at the system 

(H 1 (h)-E 1 )u = 0(h ao ), (H 2 (h)-E 2 )u = 0(h°°), (26) 

where E\ and E 2 are real numbers. Here we shall assume that Hi and H 2 
are formally self-adjoint. If we are only interested in studying the spectrum 
in a window of size 0(h) around the orig in, we can let Ei — fiEi and there 
is nothing more to be done: Theorem |2.l| holds uniformly for (ei,e 2 ) varying 
in a compact set of R 2 , so that all the results of the previous sections apply. 



However, Theorem 2A does not apply to the system (|26| ) with uniform estimates 
for E = (Ei,E 2 ) in a compact. Indeed, it would imply that H p has a unique 
value on the local level set (Hi, H 2 ) — (Ei, E 2 ) near 7^; if E is a regular value of 
F, this level set may fail to be connected and it is easy to construct a situation 
where H p has different values on each component. Actually, H p by definition is 
a function on the set of leaves of the Lagrangian foliation defined by F; and the 
following diagram is in general non-commutative: 

O— F -+U (27) 



Instead, we need to work with the space of leaves Q, equipped with the momen- 
tum map F. 

For any E G U C (R 2 ,0), let A E = F- 1 (E). If U is a sufficiently small ball 
around some critical value, the curve C c C U of critical values of the momen- 
tum map F separates the set of regular values in U into two simply connected 
open sets U + and U~ . Using H q defined in section |l,6^ we take the following 
convention : U ± := {±H q > 0}. Let D ± := U± U C c . Let N + and N~ be the 
sets of connected components of the open sets F~ 1 (U + ) and F~ l (U~) respec- 
tively. In each of U + and U~ , the levels sets of F have a unique topological 
type, namely they are unions of a finite number N of Liouville tori: for any 
Ee U±, F- l (E) =: U± eJ v±T fc ±(£). Of course, A>± = |JV±|. 

Proposition 2.8 A smooth function K commuting with Hi and H 2 in Q is 

characterised by the data of \N + \ functions /&+ € C°°(D + ) and \N~\ functions 
f k - G C°°(D-) such that 

1. For all fc± G iV*, K^ k ± = f k ± o F^ k ± ; 

2. For all k + G N + and k~ G N~ , the function egual to f k + in D + and to 
f k - in D~ is smooth on U . 



Definition 2.3 The space of smooth functions in SI commuting with Hi and 
H 2 will be denoted by C°° F (Q). The space of leaves together with the smooth 
structure described in the above proposition will be called the Reeb graph of F. 

Proof, a) Given K, the condition 1.) uniquely defines the functions f k ±. 
Their smoothness is given by the Darboux-Caratheodory theorem that states 
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Figure 7: The Reeb graph of F 



that near every non-singular point of A^, there is a canonical chart (x,y,^rj) 
in which H i = £ and H2 = r). The fact that such non-singular point exist even 
on a singular leaf Ae c , E c € C c , shows that the smoothness extends to D ± . The 
same argument shows that the condition 2.) holds whenever k + Ok 7^ 0. Then 
condition 2.) without this restriction holds because Ae c is connected. 

b) Conversely, the data of all the functions defines a unique function 
K . The smoothness of K outside of the critical points comes from the same 



argument as above; its smoothness at critical points comes from Theorem 1.5 



and the following lemma 2.£ (proved as in 123]). □ 



Lemma 2.9 Each function K commuting with £ and yrj near (x, 0, 0, 0) is char- 
acterised by two functions /+, /_ € C°°(IR 2 ,0) satisfying 

/-(£,*)-/+(£,*) = 0(*°°)> (28) 

locally uniformly in £, such that 



f-(€,m) «/y<o 



Remark 2.9. The distinction between the functions / + and /_ is of course 
irrelevant when all the data is analytic. Neither has it any impact for semi- 
excited regions, ie. for E = 0(/i 7 ), < 7 < 1 (see |2(| chapter 5]). In these 
cases, no further modification of the results of the previous sections are required 
(the diagram (27) is always commutative), and the following section p. 4.1 be- 



comes rather straightforward. We have laid down all the details to cope with 
the C°° case, which makes the statements and proofs more technical. On the 



other hand, the statements in section 2.4.2 are non-trivial even in the analytic 



case, and represent some of the most crucial results of this article. A 



2.4.1 The "global" quantum number 

In this section, the issue is to generalise the "global" quantum number e±(h) 
of Theorem 2.4. In the smooth, non-analytic category, this leads to a subtle 
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repartition property for the semi-classical spectrum. 

For any EE U\C C , denote by (£, A#) the sheaf of germs of microlocal so- 



lutions of ( |26[ ) on Ae- We know from the regular theory (see p6|) that (£, Ae) 
is just a flat bundle characterised by its holonomy, which was called in [^6| the 
Bohr-Sommerfeld cocycle, and denoted by Xh G H 1 (Ab, M/27rZ). When E is re- 
stricted to any compact subset K C U ± , and A^ is restricted to some connected 
component k , h\h has a uniform asymptotic expansion in C°°(K) This is 
no longer true on D ± . However, the following statement holds : 

Theorem 2.10 The function that assigns to a leaf T k ±(E) the integral 



h 

2^ 



A/i, 

Ft 



where 5 is any S 1 -orbit in T k ±{E), defines an element ti(h) G C°° (ft) that 
admits an asymptotic expansion of the form : 

oo 

t 1 (h) = Y J h'^ ) EC°° F (Q)lhl 

1=0 

with = Hp, and e± = J? T\ — J g r 2 + ^t(<J)/4, where Vj is the sub- 
principal symbol of Hj, and a, b in C°° F (f2) are the functions defined by X p = 
aXi + bX 2 , and fi is the Maslov cocycle of T k ±(E). 

Remark 2.10. Fix k € N^ 1 and realise (ei)ffc± as a smooth function E —> 
Zi(E;h) in C co (D ± ). Now let E = (a ,b ) + h{e,\,e2), where ej vary in a 
compact (recall that (a ,b ) is a critical value of F); then Z\(E;h) admits an 
asymptotic expansion in powers of h whose coefficients are smooth functions of 
(ei, e 2 )- But these coefficients are made out of the Taylor series of E — > ei(E; h), 
and therefore, in view of the definition of C°° F , they do not depend on the 
component k . We obtain this way an element t i{h ) £ C°° \h\ that is nothing 



else but the global quantum number of Theorem 2.4. 

Notice that, in order to compute the principal term in e\{K), it is not inter- 
esting to use the formula e\{h) = ti(E;h), since it would involve the deriva- 
tive of ej° . Instead, apply the formula of the theorem above to a system 
whose principal symbol is independent of E, and whose sub-principal symbol is 
(ri — ei,r2 — e-i). We obtain this way the claim (|lO|)-(|ll|) in Theorem 2T. A 



Proof of the theorem. The fact that ti(h) S C°° F (tt) is obvious form the 
construction. To prove the existence of the claimed asymptotic expansion, it 
suffices to microlocalize near a critical circle jj. 



Using Lemma |2.9| in the coordinates of Theorem 1.13, one sees that the 



functions / + and /_ of Lemma 2.9 are the same if the degree d = 2. In this 
case, Theorem [2.l| generalises to 

U-^H, - E 1 ,H 2 - E 2 )U = M.{Q X - t x ,Q 2 - t 2 ) + 0(h°°), 

where A4, ti, and e 2 depend smoothly on E, which gives the result. 

The case d = 4 is more intricate, and follows from Proposition 2.12 below. 

□ 



We shall need the following slightly weaker version of Theorem 2.1 
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Proposition 2.11 There exists an elliptic Fourier integral operator U(h) asso- 
ciated to the canonical transformation tp of Theorem 1.1^ such that, microlocally 
near jj : 



U-'HjU 



where Kj £ £^(7^). 



Proof. The same proof scheme as that of Theorem 2.1 applies, using Lemma 
2.9 instead of Lemma 2.2. □ 



Proposition 2.12 Let jj be a critical circle of degree 4- For each E close to 
zero, the set of microlocal solutions of ( p^ ) on a small neighbourhood of any 
point of 7j is a free <Ch~module of rank 2. In the coordinates of Theorem l.l\ 
it has a basis of the form 



Hi 



±^ e i C fx/h [ l±y>0 __^fin\y\/h 



where = tj(E\ h) admits an asymptotic expansion in non-negative powers of 
h whose coefficients are smooth functions of E. For this basis, the system $2(\) 
is solved locally uniformly with respect to (E\,E2) near (0,0). Moreover, the 
functions [E\, E2) — > — tj are flat on the set C c of critical points of F. 

The proof of this proposition relies on the following lemma : 
Lemma 2.13 Let p £ C£°(T*S 1 x T*R) be a Hamiltonian satisfying 

P(x,£,y,v) = 0fory^0 

(i.e. p—ly^Qp). Then 

Op™(p) = OpY( P ) o 1 v<0 - l y ^.Op^{p) (mod 0{h™)). 

Of course the symmetric result (with respect to y = 0) holds. 
Proof. Recall that Weyl quantisation of p is defined by : 

v(x,y) = Op^(p)u(x,y) = 



We prove the first estimate by showing that \\v\\ = 0(h°°) whenever u = l y ^ou. 
For this, we consider the two regions \y\ ^ h Sl and \y\ < h 52 , with < 5% ^ b~\ < 
1. If I y I ^ h S2 , then only the domain \^^-\ < h S2 contributes to the integral; 
and the result follows from the flatness of p with respect to its third variable : 
for all N £ N, \v\ = 0(h N52 ). Let us now look at the region \y\ > h Sl . Since 
v(-, -,y 0, •) = 0, one can assume that y ^ —h Sl , which implies \y — y'\ h Sl . 
Now the usual trick applies : a repeated integration by parts with respect to 
the operator w^z^n ( or standard estimates for the Fourier transform) gives 

\v\ = 0(h N ^- Sl ^) for any integer N. 
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The same method can be applied to show that ||lj,^cH| = 0(h°°) whenever 
u = l y ^o u , thus proving the second estimate. □ 



Proof of Proposition 2.12. The fact that the set of solutions is a free module 
of rank 2 is due, for E $ C c , to the regular theory (the local Lagrangian manifold 
has two connected components, on each of which the set of solutions is a free 



module of rank 1), and, for E G C c , to Proposition 2.5. 

We prove the rest of the proposition for u%; the same argument applies to 



u E . First apply proposition 2.11 to assume in what follows that Hj € *£-h{lj)- 
Since yrj is a quadratic function, every element of if € £/ l (7 J ) can be written 
K = Opf{K h ), with K h ~ J2i>o h £ KW and g £ d ( 7i ). Because of Lemma 
2.9, each is defined by two functions /+. Let 



± 



and Rh = Kh — Fh- Let us prove now that there is a unique symbol gh{£,,t) ~ 
E/^o^V^Cfi*) such that 

OpF(F h )=g h (Q 1 ,Q 2 )+0(h 00 ), 

(Qj is defined in Eq.(|J)). Indeed, g^ is necessarily equal to f+; therefore, 

Opl v (F h )=gM(Q 1 ,Q 2 ) + hS, 

where S € ^-h{lj) an d is of order 0. Then S is similarly decomposed - and the 
claim is proved by induction - provided we show that its Weyl symbol is, as for 
Fh, a function of (£, yrj). This is achieved by applying Lemma and remarking 
that Op\ v (Fh), as well as g^ (Qi, Q 2 ), commute with the involution y — ► —y 
(and thus their Weyl symbols are invariant under (y, 77) — > {—y, —rj))- Summing 
up, we have proved so far that any operator K G €h{"fj) can be written : 

K = g h (Q 1 ,Q 2 ) + Op h v (R h ), 
where all the coefficients in the expansion of Rh are smooth functions verifying 



the hypothesis of Lemma 2.13 . 

Applying this to Hj, we obtain the existence of two symbols g±.h and g 2> h 
such that 

H jU + = gj , h (zi,4)u+ + 0(h°°). 
The independence of Hi and H 2 implies that the principal term : 

(tt)^(g^(tt),g { 2 0) (U)) 

is local diffeomorphism; therefore the symbol (gi,/i, g 2 ,h) is invertible, and the 
proposition is proved with 

{tt,ti)~{ g i, h ,g 2 , h )- x {E u E 2 ). 

□ 
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Corollary 2.14 Fix fc* € AT* and realise (ti)^ k ± as a smooth function E — > 
ti(E;h) in C oa (D ± ). Then the condition 

ti{E;h) E KL + 0(h°°) (29) 

is necessary and sufficient for the existence of a uniform solution of ^2(\ ) mi- 
crolocalised in a neighbourhood (in k ) of any S 1 -orbit in k . 



Definition 2.4 For any natural integers £ + and £~ , we call an {£ + , £~)- curve 

the union of £ + smooth curves in D + and £~ smooth curves in D~ that are 
transversal to C c and infinitely tangent to each other on C c . 

These curves are just the image by F of a level set of a smooth function K on 
the Reeb grap h of F, if Xr ^ on the critical leaves A# c (this is a consequence 
of Proposition 2.8). This holds for instance for K = H p (Fig. ||). 



Corollary 2.15 There exists a fixed neighbourhood U in ffi . = {(E±, E2)} of 
any critical value of F in which the joint spectrum of Hi and tl% is distributed 
(modulo 0(h°°)) on the union of (\N + \, \N~ \)-curves L n (h) (n £ 2,) defined as 
the image by F of the level sets ti(h) = hn. The principal part of these curves 
is thus given by the level sets of H p . 




Ei 



Figure 8: Level sets of i? 3 



Proof. Fix k E N . Then Proposition 2.8 says that the restriction (ei(h)) tT ± 

I 1 I k 

is equal to /j,oF-± for some smooth function f^ ~ J27Lo ft admitting an asymp- 



totic expansion in C°° (K , 0) [ft.] . Since /oof -± = (H p )^±, the hypothesis (|6|) 
of Theorem 1.6 implies that the foliation f^ — const is transversal to C C1 and we 
can define the projection vj^ : K 2 — > C c such that vjh{f = const) is a point. The 
pre-image of {A = const} DD ± by F tT ± is a leaf of the foliation {ei(ft) = const} 

in k . The value of H p on C c can be taken as a coordinate on C c , and via this 
identification, it is natural to view ei(ft) as a function with values in C c . By 
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Corollary 2.14, any microlocal eigenfunction of {H.\,Hz) microlocalised in k 
defines a joint eigenvalue that belongs to w^ x (JiL + 0(h°°)). 




(30) 



□ 



2.4.2 Regularization of Xh(E) 

This section contains some of the most central results of this article (Theorems 



2.19 and 2.20). They give a new interpretation of the holonomy of section 2.3 



and provide for the proofs of the various formulae claimed before. 

To each kr S A* and E c £ C c we associate the subset T k ± (E c ) = k D Ae c 
which is the "limit" of the torus T k ±(E) as E — > E c . Suppose we are given 
a continuous family of piecewise differentiable loops (Se)e^d ± on T k ±(E) that 
are everywhere transversal to the S^-orbits (such a family can be constructed 



using for instance the normal form of Theorem 1.13 ), and assume that they are 
oriented by the flow of H q . For non-singular values of E, Se together with an 
S' 1 -orbit form a basis of Hi (T fe ± (£")), and E — ► f s Xh is a smooth function. 
To complete the result of Theorem [2.10] , it is natural to investigate here the 
behaviour of that function as E approaches a critical value. 

To each E c S C c corresponds a real number x via the diffeomorphism (|30|); 
T k ±(E c ) is an S' 1 -invariant subset of H~ (x), and hence can be reduced to a 
cycle G k ±(E c ) of the graph G(E C ) :— G(A Ec ) in the reduced manifold W{E C ). 
The second goal of this section is to show how the asymptotic behaviour of the 



function E — > J s Xh is related to the holonomy hoi of section 2.3. Actually we 
shall prove that J s Xh diverges as E approaches a critical value; but there is a 
universal way of regularizing this divergence. The regularized value is precisely 
Yio\(8e) ~ modulo some Maslov corrections in presence of vertices of degree 2 
-, where Se is the projection of Se onto the reduced manifold and is actually 
equal to G k ± (E c ). 

Unfortunately, the set of cycles G k ± (E c ) doesn't necessarily generate the 
group Hi (G(E C ), Z) - see eg. Fig. |9|- (but it does indeed if G{E C ) is planar). 
So there's a little bit more to it than just taking the limit of regular cycles. 

To compute hoi of all possible cycles of G(E C ), we have to replace the 
"natural" object Se with a local path near a critical circle. In order to give 
some sense to the expression J s Xh, where Se is not closed, we could abstractly 
choose a smooth family of closed 1-forms on A^, E G U , whose cohomology 
class is [Ah(-E)] - which is always possible as H 1 (Ab) — ► U is a trivial bundle 
if U is small enough. However, this does not allow us to have a local control of 
the divergence of the holonomy. Instead, we follow [^7j and interpret A;,, as the 
phase of the multiplicative Cech holonomy of the sheaf (£,Ae), as follows: 

Definition 2.5 Near each Tj, we let Aj(E) and Bj(E) in C co (D ± ,k ± \ Tj) be 
families of points such that for a critical value E c , A(E C ) and B(E C ) lie respec- 
tively in the local stable or unstable manifold. We endow a small neighbourhood 
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Figure 9: A graph on a torus. The cycle in gray cannot be obtained as a 
combination of boundary faces of the graph. 



of Aj (or Bj) with a standard WKB microfocal solution ua- (resp. ub ) whose 
phase admits an asymptotic expansion in h~ 1 C°°{D ± )\h\. 

Then the integrals J A 3 ' \h are defined as the phase of the Cech holonomy of 
(£, K-e) for paths joining Aj and By with the sections ua^ and Ub , fixed. 

In other words, if the path 5 between Aj and Bji is covered by open sets 
flo, ■ ■ . each of which being endowed with a microfocal solution u a with 
uq = ua and ug = ub ., , then 



J A/j = -i In (co,iCi !2 • • • Ci-\,i) , 



where Cjj- is the transition constant Ui = CijUj on Qi PI Qj. 

Note that if another admissible choice for the local sections ua/b 1S made, 
then the holonomy is modified by an additive term admitting an asymptotic 
expansion in h~ 1 C°° (D ± )lhJ. Therefore the singular behaviour of the holonomy 
at a critical value is fully preserved. Note also that this additive term is necessary 
a Cech coboundary, and hence has no influence on the value of the holonomy 
along a closed loop. 

Definition 2.6 In what follow, {Se)e£D ± designates a continuous family of 
paths in T k ± (E) such that: 

• for each E € U , Se is smooth; 

• either for all E c £ C c , Se c does not meet the critical set T (then (6e) is 
called regular) or for each E c , 5e c meets uniquely a unique critical circle 
jj (in which case (Se) is called local); 

• the end points Aj(E) and Bji(E) are one of those defined in Definition 



2.5, and we will write Se = [Aj(E), Bji (E)\; 



5e is always transversal to the S 1 -orbits. 



Here again, the normal form of Theorem [1 . 1 3] proves the existence of such a 
family of paths near any critical circle jj . 

The goal of this section is finally to investigate the behaviour of the function 
E — > J Se A/j as E tends to a critical value, and to relate it to the holonomy 
of the sheaf (£, G(E C )). The previous case where Se was a loop can always 
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be recovered by composing paths of the type of Definition |2.6| . Moreover, the 
regular theory implies that the so-called "local" paths can indeed be restricted 
to paths that are local in small neighbourhoods of the critical circles, since the 
following proposition holds: 

Proposition 2.16 If Se = [Bj(E) > Aj/(E)] is a regular family of paths (in the 



sense of Definition 2.6), then E — ► f s Ah is smooth in D ± and admits an 
asymptotic expansion in h~ 1 C°°(D ± )lhJ. The first terms of this expansion are 
the following: 

[ \ h = * Bl {B j )-* A {A r ) + \ [ a+ [ K + ix{5 E )~+0{h), (31) 

JS E n JSe JS e 1 

where $>Bj (resp. ^A jt ) is the phase of the principal symbol (viewed as a section 
of the Keller-Maslov bundle over the Lagrangian manifold A E \^j - see eg. Jlq / 
or °f the fixed solution ub } (resp. ua., )■ 



To study the neighbourhood of a critical circle, we shall use Theorem 1.5 
to generalise the semi-classical invariant £2,j{h) of equation ( |l7| ) in a better 



way than Proposition 2.12 would do. That theorem still holds if Hi and H 2 
are replaced by Hi — E\ and Hi — E2, for a parameter E — (Ei, E2) varying 
near (a ,b ). If wc fix a critical circle 7^ and m <E Jj, the theorem yields a 
canonical change of coordinates (x,y,£,,rf) near m, depending smoothly on E, 
and a function <j># G C°°(R 2 , 0) depending smoothly on E, such that 

H l -E l = e, H 2 -E 2 = <f> E (£, yrj). (32) 

(We have still 82^ e (0,0) > 0.) This leads to yet another semi-classical normal 
form : 

Proposition 2.17 Let 7j be a critical circle and m € Jj- There exists a 
microlocally unitary Fourier integral operator U(h) associated to the canoni- 
cal coordinates (x, y, £, rf), elliptic pseudo-differential operators Mi(h), M 2 (h) 
commuting (modulo 0(h°°)) with Q\ and Q2, and a real-valued function of 

h (independent on m, U, Mi): z 2 = ^2 tj (E;h) ~ Efco^j^^* such that ' 
microlocally near m : 

U-^m-E^U = Qi + 0(h°°) 

U-\H 2 -E 2 )U = M 1 Q 1 +M 2 .(Q 2 -e 2 ) + 0(h QC ). 

Mi, U and t 2 depend smoothly on E. 

• t 2 °\{E) is equal to the value oj 'yrj on Kb- In particular, e^ ] > if E G U + 
and 4°j <QifEeU~; 



• IfE = E c G C c 



where Ti is the sub-principal symbol of Hi, A = X(E C ) is the unique real 
number such that H 2 — \H\ is critical on Tj fl (see Lemma and 
\H.y;(H 2 )\ is the absolute value of the determinant of the transversal Hes- 
sian of H 2 . Note also that this denominator is equal to 8 2 <S> E (0, 0), and 
A = 9i$ E (0,0). 
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Remark 2.11. (See Remark 2.10). If E is restricted to a domain of the form 
E = (a ,b ) + h(ei,e2), where (a a ,b ) € C c , then Theorem [O] applies with e\ 
and ei as parameters, and yields an invariant C2j, which can be recovered from 
62 j by the following formula : 

£2,i(ei,e 2 ) = ^t 2 ,j((a ,b ) + h(e 1 ,e 2 j) + 0(h°°) 7 
or merely by viewing — (ei, e^) as a correction of the subprincipal symbols and 



applying the formulae of the Proposition. This proves the claim (12) of Theorem 



2T| A 

Using this proposition, we let (3 = ^+ie 2 /h and be the /i-dependent functions 
in C°°{U) defined by 

C+ := - 1 L=T(l3)e^ h e- i ^=e- i i£ j (t2,j/h); (33) 
q := - 1 L=T(f3)e^ nh e 1 ^ ^e l Te-^ /h £ 3 {e2,j/h). (34) 



/2nh 

(£j was defined in (^).) Next lemma, which directly follows from Stirling's 
formula, will be crucial for our analysis. 

Lemma 2.18 For any E e U ± , 

iln(f = i(egln|eg|-eS) + eSln|eS| T J + 0,(/ l ). 

Theorem 2.19 Fix a component kr 11 , and let 6e = 



Aj(E),Bj(E)] be a local 
path near the critical component Tj (see Definition 2.t ). Assume moreover that 
Se is oriented according to the flow of H q (otherwise just take the opposite 
of the holonomy !). Then there exists an h-dependent M./2n'Z-valued function 
gs(h) : E — ► gs E (h) E C oa (D ± ) admitting a uniform asymptotic expansion of 
the form 

oo 

gs(E; h)~J2 98 } (E)h e , gf G C°"{D±), 
i=-\ 

such that 

V-EeC/*, g 5 {E;h)=[ X h - iln((f (E)) (mod 2vrZ). (35) 
The principal terms of gs(h) are given by the following formulce, for E G U ± : 
9^\E)=j a+ (egln|42|-eg)+^(^-^-. 1) (5,); (36) 



9 { s°\e) =f ia K*+ + (t \ + 43 ^ i43i) + <W - *%m 

(37) 

where ^>^ j = + &a is the phase of the principal symbol of the fixed 

solution UAj (and similarly for UBj)- 
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Theorem 2.20 Let E c — (a ,b ) G C c . Let 5e c be a loop in Ae c oriented 
according to the flow of H q and of the form 

X _ — Xk>c rrcg doc rreg Hoc rrcg 
ffic — °\ °\ 2 2 O q O q , 

where c^ oc and 6^ are respectively "local" and "regular" paths in the sense of 
Definition \2. (\ . (The components fc* used for these paths may vary. - see eg. 
Figfy. Let 



g(E c ; h)~J2 9 W (E c )h e , S C°°(C C ) 



be defined as the sum 

g{E c ; h) := (g s u,c + g s *« + ■■■ + g s i»c + g s ^ ^ 



where g s iac is given by Theorem 2.H and gg^s d = fj le g A/j (see Proposition 2. It ) 



Then, under the hypothesis of section 2.1 



ho\{8 Ec ) := e i9(Ee;h) e ifn(N 2 -N+) + o(h°°), 



(38) 



onto the graph G(E C ) in 



where Se c is the projection of Se c 

ifold W(E C ), n is the "global quantum number" of Theorem 2.4, and N 2 



the reduced orb- 
is the 



number of local paths through a vertex of degree 2 and defined by a component 

fee N ± . 



Proof of Theorem 2.1(\ . Fix a critical value E c S C c . Let m be the inter- 
section of the cycle 5e c with jj, and Q m an open set in which Proposition 2.17 
applies. We can assume that the paths Se, E G D all entirely lie in fi m (using 
Proposition 2.16, this will only modify J s A/j by an additive term entering in 
gs(E;h)). As before, label the local un/stable manifolds with cyclic order I lt 
7 3 , I 2 and I4. 5e c enters tt m on a local stable manifold I s , s — 3,4, and leaves 
it on a local unstable manifold /„, u = 1, 2 (it is the index preceding s in the 
cycle (1, 3, 2, 4) if k^ 1 = k + 6 N + and the index following s if fc* = k~ E 
As before, we endow a neighbourhood of each I a with the distribution e a : 



P4 



• — J -±y>0 f — |-o 



1 



1 , .„ Z P «2,jln|j/|//i 



^ 1 («: ,a, ' /fc )- 



These distributions are classical Lagrangian distributions whose phases admit an 
asymptotic expansion in j^C°°(D ){hj. Moreover, they are microlocal solutions 
of ( p6|) in Q m unifor mly for E € D ± , and hence constitute an admissible choice 
in view of Definition 2J). Note that this choice possibly implies another additive 
term entering in g$ E . In a small ball around Aj(E c ), the space of solutions has 
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dimension 1 and we must have a constant C ± (E, h) such that e u ~ C ± (E, h)e s 
Now Definition 2.5 says that, with respect to the fixed solutions e u and e s , 



= e 6 



(mod 0(/i°°)). 



The expression of the Fourier transform involved in e s shows - as in [pi - that 



C ± (E,h)=d 



(mod 0(h°°)). 



(39) 



This proves the existence of the claimed gs satisfyin g Eq. ( p5| ). The formulae 
( |36"1 ) and ( |37| now follow from Stirling's formula (Lemma 2.18 ) and from the fact 
that the asymptotic expansion for gs is uniform in D , and hence is given for 
any fixed E e U by the difference of the (non-uniform) asymptotic expansions 
of J s Xh and i m (Cf(E)). We have also used the formula ([II]) which holds 
only for non-critical E, and which comes from the definition of the bundle of 
principal symbols. 

Note that we don't see here any difference between critical circles of di- 
rect/reverse type, for the statement is purely local near a point of jj. □ 

Proof of Theorem 2.2C. To each S\ oc is associated a unique component k se , 
Si = ±, a unique critical circle 7^, and a (non- unique) local family of paths 



5f c (E) (in the sense of Definition |J) such that df c = Sf c (E t 
set of E £ D s * such that (eOrfc»* G hZ + 0(/i°°) (ie. T h l = w^[hZ + 0(h°°)} 
in the notation of the proof of Corollary |2.14j ). Then the (assumed) hypothesis 
of section 2.3 "ex = ri" says that E c S 2^. Using Proposition 2.12| one can 



Let T h l be the 

v3l x \KL- 



construct smooth families (ue) 
k Sf such that 



1,2,3,4 
Eexf e 



of solutions on a neighbourhood of jj e in 



(u E ,u 



El a Ei 



4) 



(40) 



the direct case - or 



e J 7"ui 



(41) 



- in the reverse case - form at E = E c a standard basis for the graph G(E C ) at 
the vertex gj e . 
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Since these solutions are smooth WKB solutions and hence admissible in 



the sense of Definition 2.5, we shall use them to define the local holonomies 
J^ioo/teg^j A^; since they are standard basis at E = E c , we shall in the same 

loc /rcg 

way use them to define the local "reduced" holonomies hol(<^ ). But then 
by definition of the sheaf (£, G(E C )) we have 



hol((^ oc ) = 1 and 
hol(^ cs ) = exp (ijg** 



exp 



(ig s ™z(E c ) 



(42) 



On the other hand, we know from the proof of Theorem 2.19 that for such 
a choice of microlocal solutions, we have (modulo 2tt) 



g s iac (E) = if jj t is of degree 4; 
g s iac(E) = se^n if -fj e is of degree 2. 



(43) 



Therefore, if we decompose 



hol(5 Ec ) = [] hol(,5 

£/ deg 73f =2 

we obtain by @ and (p3h : 



—loc 



x n hoi^° c )xnhoi(C), 

£/dog 7jf =4 e 



hol(d Ec ) = exp ii (5«5>° c i. E c)-H^i 

\ e/dc g7je =2 



J2 g sr {E c ) + Y.9si<E c ) 

' dog 7j =4 I 



which proves the theorem. 



□ 



Corollary 2.21 Theorem 2.2C together with formulae j\36\ ) and ^3\) finally prove 
the second point of Proposition 2.6. 



Proof. 



1. The principal action — Since t^{E c ) = 0, it is clear from (|36|) that for 
any 1= 

(E c ) = [ a + (Aj e (E c )) (B je (E c )) 

and from ( ^l] ) that for any £ = 1, . . . , q (identifying I = q + 1 with £ = 1), 
gg? (E c ) = f a + (B je (E c )) - $™ {B h+1 (E c )). 

Therefore, 



g { - 1] {E c )= I n. 

Isr 
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2. The sub-principal action and the Maslov index Fix £ = 1, 

and let 7j = 7^ . Let m = m(E c ) be the point where <^ oc meets 7,. As 



Proposition 2.17, we shall use the local canonical coordinates at m given by 
Theorem 1.5 



Recalling the notation of Proposition 2.17, (B2|) implies that 



1 



(x 1 ,x 2 )-[ d ^ E d ^ E )-(x e ,x l 



Since 82® e ^ 0, there exist a smooth function p 2 = dl d^^ +T2 (depending 
also smoothly on E) such that the sub- principal form ke is given by 

k E -{X£,X vv ) = -{r ll p 2 ). 

Note that for a critical value E = E c , (p2)\ y=n =a = — eg- The closedness of ke 
on each implies that 

{ri,yv} = {P2,£}- 

Using a local analogue of Lemma |2.3| , we can decompose (r*i, p 2 ) in the following 
way: 

(r 1>P2 ) = (0,K) - (44) 

for some smooth functions /, K where K commutes with yrj and £. Therefore 
the function 

/ := / — K In \y\ (or / + AT In \r)\ where y = 0), 
restricted to A#, satisfies d\ E f = Ke- We can now compute 

/ Ke — f(B 3 ) - f(Aj) = f{B 3 ) - f(A 3 ) - K In \y B]VA] \+Kln \y AjVAj \ 
Since A 3 {E) and B 3 {E) are in A^, VAjVAj = eg an< ^ we nave 



f K E + eg In |eg| = f(Bj) - f(A 3 ) -Kin \y Bj VA 3 \ + (eg + K) ln|e£ 



Because of (f44|), pi — K vanish at y = r\ = 0, hence K = —eg + 0{yrf) and 

(eg + K) In |eg | tends to zero as E tends to a critical value E c . Since g^} c is 

smooth at E c , the formula ( |37| ) implies that /i(<^ oc (i?)) is continuous at E c and 
hence constant; let us denote it by p(S\ oc (E c )). 

Suppose now that E = E c and let a, b be points on 5\ oc (E c ) located respec- 
tively in [j4j,m] and [to, Bj]. Then 

/(^)-/(^)-A-ln|y fl ^| = 
= lim (/(a) - f(A 3 ) + f(B 3 ) - f(b) - Kin \ yBj rj Aj \) . 

a,b— >m \ / 



2 SEMI-CLASSICS 



39 



The term in the limit is equal to 



[A,, a] J[b,Bj] 



Therefore, by Definition 2.2, the limit is equal to J 5 i oc ke c , and 



yields: 



Then as before, if we sum up all the contributions from regular and local paths, 
we obtain 



,(0) 



Using the definition 2.2 of the regularized Maslov cocycle, we finally obtain 



g^(E c ) + £(JV 3 " _ N +) = kEc + ~^5 E f-, 



and equation (B8h concludes the proof. 



□ 



3 EXAMPLES 



40 



3 Examples 

We propose in this section several examples for which our theory applies. Many 
other could probably be found; the ones presented here are interesting by their 
simplicity and yet by their rich structure and behaviour. 



The recipe 

Let us recall here briefly the recipe for obtaining the semi-classical quantisation 
rules. The first thing to do is to locate the critical value of transversal hyperbolic 
type in the image of the momentum map F — (Hi,Pi2). Then choose one of 
these points o and describe the singular level set A G = i 7 '~ 1 (o), in order to have: 
a) the graph G, b) a formula for the vector fields X\ and X 2 on A Q - and if it is 
not one of these, for the periodic vector field X p . 

Compute the semi-classical invariants (action integral, Maslov index) for a 
periodic cycle - this implies only reg ular tools - in order to derive the first 



quantisation condition of Theorem 2A up to 0(h), and fix the quantum number 
n. 



From the graph, apply Theorem 2.7 to obtain the second quantisation rule in 



the form of a determinantal equation. It remains to compute the holonomy hoi 



up to O(h), which involves the singular semi-classical invariants of Definition 2.2 . 
The fulfilment of these quantisations rules determine the spectrum up to an 
error of order 0(h 2 ) in a window of size O(h) (in fact, it is easy to determine 
the smooth dependence of the semi-classical invariants in o, and the spectral 
window can be extended to a rectangular domain of size O(l) along the curve 
of critical values, and of size O(h) in the transversal direction). 



Notation— The reader must be warned that the symbol e (italic) is used as 
a subprincipal spectral parameter (as in " E — he" ) , while the exponential is 
denoted by e a = exp(a). 

3.1 Laplacians on Ellipsoids 

Let us consider the ellipsoid in the Euclidian space R 3 defined by : 

9 ? 2 

E = {^ + ^ + ^ = l} 
a{ a\ a\ 

with < a\ < a 2 < a 3 . The geodesic flow on E has been discovered to be 
integrable by Jacobi in 1838 using Abelian integrals. For a recent presentation, 
one can read §|, §, @ or ||. 

3.1.1 Classics 

Let us denote by P,Q,P' = —P,Q' = —Q the four umbilics of E which are 
located on the ellipse {x2 — 0}. If 

X\ = v / «i(a2 - a 1 )(a 3 - Oi) and X 3 = ^a 3 (a 3 - a 2 )(a 3 - 0%) , 



we have 



P=(X U 0,X 3 ), Q = (-X 1 ,0,X 3 ) 
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We will consider the (unique up to global dilatation) conformal representation 
$ of E + = E n {x 2 > 0} on a rectangle R =]0, Ti[x]0, T 2 [ such that the four 
umbilics are going on the four vertices of R, according to figure [l4|. Using such 
coordinates (x,y) <E R, we get (sec Q (vol. 2 p. 308 and vol. 3 p. 13) or ppf) 
the following expression for the metric of E : 

ds 2 = (a 2 {x)+b 2 (y))(dx 2 + dy 2 ) (45) 

where a, b are given in terms of hypcrclliptic integrals and extends to smooth 
functions on K which satisfy: a is > on ]0, Ti[, vanishes exactly at the points 
kTi, k E It, and is odd with respect to T\1 and b satisfies the same properties 
with respect to T 2 . Moreover a'(0) = b'(0) > 0. Let us denote by T the lattice 
Til © T2Z. Then ds 2 extends into a smooth metric on R 2 \ T, which is L- 
periodic. Let us consider the torus T = R 2 /2L. Then the map a : T — > T 
defined by T(z) = —z defines an isometric involution of T with four fixed points 
and we get a natural identification of E with T/a as a 2-sheeted branched 
covering II of T over E with automorphism a. The metric ds 2 admits conical 
singular points of total angle 47r at the umbilics which makes the metric on E 
smooth. More precisely, it follows from the formulae of |2(| that there exists 
an analytic function G defined near 0, with G(0) = 0, G'(0) > such that 
near (0,0) we have ds 2 — (G(x 2 ) — G(—y 2 ))(dx 2 + dy 2 ). It is rather easy to 
check that there exists an analytic function A(u, v) with A(Q, 0) > such that 
G{x 2 ) - G(-y 2 ) = {x 2 + y 2 )A{x 2 - y 2 , 2xy) and, if locally Z = TL(z) = z 2 , we 
have : ds 2 = Ii*(AA{X, Y)(dX 2 + dY 2 )). 

We will use the fundamental domain D = [0, 2T±] x [0, T2]. E can be recov- 
ered from this rectangle by gluing edges as indicated on the figure ^4| 

If a, b were non vanishing, ds 2 would be called a Liouville metric on T. Our 
case corresponds to a degenerate Liouville metric on the sphere. It is well known 
that Liouville metrics are integrable. Let us denote by 

Hi= e+v 2 



the geodesic flow and by 



H 



2 



a 2 (x)+b 2 (y) 
b 2 (y)e -a 2 (x) V 2 



a 2 (x)+b 2 (y) 

The manifold Le,f = {Hi = E, H2 = F} is given by: 

Le,f = {i 2 = F + a 2 (x)E, V 2 = b 2 {y)E - F} 

which is obviously Lagrangian. 

We are interested in the singular value o — (E = 1, F = 0) of the moment 
map (Hi,H 2 ) and the corresponding A D . Geodesies passing through P (resp. 
Q) contain also P' (resp. Q') and vice-versa. A D is the set of unit covectors 
corresponding to geodesies passing through P or Q. We have A G = \JL± t ± 
where L± ± = {£ = ±a(x), rj = ±b(y)}. In particular Lp = L + + U L_,_ and 
Lq = L+,_ U £— + are smooth Lagrangian tori whose intersection is 7+ U 7_. 
Here 7+ (resp. 7_) is the lift of the ellipse x 2 = with orientation (P, Q, P' , Q') 
(resp. opposite). Lp (resp. Lq) is the set of unit covectors of geodesies of E 
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containing P and P' (resp. Q and Q'). Lp (resp. Lq) is the stable manifold of 
7_ (resp. 7+) and the unstable manifold of 7+ (resp. 7-). 

The associated graph G is the union of 2 circles corresponding to Lp and 
Lq intersecting at 2 points corresponding to j± . 




Figure 12: associated graph for E 




3.1.2 Quantum 

We now introduce the quantum Hamiltonian Hi = /i 2 A^ which is given in the 
coordinates (x, y) by: 

h 2 

Hl = ~ a 2 + b 2 (d* + dy) 
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Figure 14: fundamental domain 



and the operator 



Hi = -^—Ab 2 d 2 x -a 2 d 2 



b 2 



It is possible to check directly that [Hi, Hz] = outside of the lattice T which 
is the set of singular points. We need a stronger form of commutation, namely 
if 

Hup = Xtp , (46) 

we need to prove that H 2 tp is in the domain of Hi and H\{H2p) = XH 2 p , so 
that Hi and Hi have a commun eigenbasis. 

For that, we need a caracterisation of the functions <p = (p o II where tp is an 
eigenfunction of the Riemannian Laplacian on E. 

Lemma 3.1 tp : T — ► C is of the form tp = tp o II with Ae<P — ^<P if an d only 
if ip is a invariant and satifies 

- = \(a 2 + b 2 )(p . (47) 

Proof. Starting from tp an eigenfunction of A^ and using smoothness of II 
gives the trivial direction. For the other, using the fact that tp o a — tp, we get 
a bounded tp with tp = tp o II and (A^ — \)tp = T where T is supported insidcs 
the (finite) set of umbilics. We deduce that T is or that tp (and hence tp) is 
unbounded. □ 

It is easy to check that Hitp is smooth using equation (Q). Then Hi'p sat- 
isfies equation (|47| ) outside T and hence everywhere on T and we get a commun 
eigenbasis for Hi and H 2 - 

We prefer to rewrite the eigenvectors equations 

Hitp = Xtp, H 2 tp = [itp 
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in the following simpler way: 

P(p := h 2 -^ + (a 2 {x)\ + n)ip = Q 

Q£ ■.= h 2 -^ + (b 2 (y)X-fx)<p = ■ 

We are interested in solutions of this system which are a invariant. If we denote 
by a%(x, y) — (—x, y) and <72(x, y) = (x, —y), we get a = o\ o <j 2 and because <7j 
commutes with P and Q we are reduced to find solutions of the form (p(x, y) = 
f(x)g(y) with / a 2Ti-periodic solution of Pf = and g a 2T2-periodic solution 
of Qf = 0. We ask moreover that / and g are both even or both odd. We 
assume A = 1 which corresponds to quantize h and fi = eh. The associated 
fiber of the momentum map is then A Q . 

This way we are reduced to 2 one dimensional problems and because P and 
Q are semi-classical stationnary Schrddinger operators with potentials —a 2 and 
— b 2 , we are reduced to the computations of p. 489-490] for periodic double 
wells. 
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3.2 1 : 2-resonance 

3.2.1 Birkhoff normal forms 

Consider a Haxniltonian H : T*R 2 —> M with a non-degenerate minimum at 
the origin. We can assume using a symplectic linear change that H(zx,Z2) = 
K 2 {z) + 0(|z| 3 ) with 

K 2 (z) = uji\zi\ 2 + uj 2 \z 2 \ 2 , 

and u>j > 0. Here Zj = Xj + where (x\, x 2 , £i, £2) are canonical coordinates 
for T*R 2 . We will say that the quadratic part is resonant if u>i/u 2 is a rational 
number. It is possible to derive a Birkhoff normal form H of the following form 

H = K 2 + R+0(\z\)°° 

with R — 0(|z| 3 ) and {K 2 ,R} — 0. The same result is true on the quantum 
level (see chapter 5 of pq ) with commuting operators K 2 and R . If we are 
able to analyse the joint spectrum of the operators K 2 , R we can deduce some 
sharp results for eigenstates in the energy domain E = 0(h a ) with a > 0. In 
the case of the 1 : 1 resonance - ie. u>i = lo 2 - the flow of K 2 induces a free 
circle action on the energy hypersurface K 2 = const and the reduced space is 
smooth. Then, via the use of Toeplitz operators, the problem is fully reduced 
to a 1-dimensional one. This is no longer the case for the 1 : 2 resonance, where 
the reduced phase space has a conical singularity. For this simple example, we 
will show that our analysis applies. Another application would be the near 1 : 2 
resonance with R — e\zi\ 2 + R . 

3.2.2 1 : 2 resonance 

We consider the following Poisson commuting Hamiltonians on T*M. 2 : 

Hi = ^\ Zl \ 2 + \z 2 \ 2 , H 2 = (x\ - g)x 2 + 2x^2 = 5R( Zl 2 z 2 ) (48) 

with Zj = Xj j = 1,2. The image of the momentum map F — {Hi, H 2 ) is 

F(T*R 2 ) = {(X, Y)\ 16X 3 ^27Y 2 }. 

The singular values consist of the boundary (which corresponds to transversal 
elliptic points, except for the origin which is degenerate) and the half line C c = 
{(X, 0), X > 0}, whose points are transversally hyperbolic (see Fig.|l5|). Here 
we shall be interested in the critical values on C c . Because of the homogeneity 
of Hj, it is sufficient to consider the point o = (1,0). 

The corresponding commuting quantum Hamiltonians are: 
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Figure 15: Image of the momentum map for the 1 : 2 resonance 



3.2.3 Classical description 

Here we are interested in the singular Lagrangian leaf A Q defined by Hi = 
1, H2 = 0. The singular part of A Q is the closed trajectory 70 = {z\ = 0}n{|z2| = 
1}. From its defining equations, it is easy to find a parameterisation that shows 
that A„ is a Lagrangian immersion of a Klein bottle K with 70 as a double loop: 

$:K9(^)h (V2e 10 sin ip, -ie 2W cos ^)eA„, (51) 

where K is the quotient of T 2 e <. = K 2 /Z 2 by the equivalence relation 

A fundamental domain D is given by D = {(9,<p) | 0^9^ tt, — tt ^ 
tp ^ 7r}. The singular line 70 corresponds to {p = 0} U {p = ±ir} and we 



have there the identifications <S>(6, 0) = &(9 + n/2,±ir). The graph G = 




corresponding to A Q has just one vertex 70 and one edge. 

Remark 3.1. Although we don't really need it, it can be helpful to have 
a representation of the reduced phase space W = H^~ 1 (l)/S 1 - where the S 1 - 
action is the flow of the harmonic oscillator H±. Using a priori argument, one 
can show that W is a 2-sphere with a conical singular point; however, one can 
find an explicit equation for W. The algebra of ^-invariant polynomials - that 
is, those that commute with Hi - is generated by 

7Tl = \zi\ 2 , 7T 2 = |Z 2 | 2 , 7T3 = ^(^2), 7T 4 = j(zfz 2 ) , 

which are subject to the relation tt 2 + ix\ — tt\tt2- One can show (see the 
book ||) that this relation together with 7Ti ^ and 712 ^ define the orbit 
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Figure 17: Parameters set for A Q 
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space T*^ / S 1 in terms of the variables Wj. The energy level set is the section 
{tti 0} n {tt2 ^ 0} of the 3-dimensional hyperplane tti + 2-k 2 = 1. Therefore, 
W is defined in the space R 3 = (tti, 7T4, ^3) by the equation 

tt| =7Ti(l-7ri)/2-7r|, with tti e [0,1]. 

PF is a surface of revolution around the 7Ti-axis, homeomorphic to a 2-sphere, 
with a conical singularity at the origin. Note that H 2 = ^3 so that the restriction 
(-^2)^ is just the height function (see Fig. |l8| ) and is a Morse function on 
W \ {0}. The manifold A D reduces to the singular equator 7r 3 = 0. A 




Figure 18: The singular reduced phase space W. Here the function (H 2 )\w is 
equal to the 773 coordinate. 



3.2.4 Semi-classical computations 

We consider the solutions of the system 

(Hi - l)u = 0, (H 2 - eh)u = 0, (52) 

for bounded e. The microsupport of the solutions is A D . 

If we denote by Xj the Hamiltonian vector fields of Hj we get on A Q : 

^ = -|,^ 2 = -2si n ^, (53) 

which leads to the following sub-principal form k for F — (Hi , H 2 — eh) : 

k = -^dip . (54) 
I sin if 

Note also that, since the flow of X\ is 2-7r-periodic outside 70, we have X p = X\. 
The canonical 1-form a = ^\dx\ + £, 2 dx 2 is given by: 

a = -2(sin 2 6»sin 2 <p + cos 2 29 cos 2 ip)dd + ^ (sin 20(1 + 008 20) sin2(p)d(p . 
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Finally, the caustic set C of A D is given by 

C = {cos 6 = 0} U {tan 2 6 = cos 2 <p}. 

In order to compute the quantisation rules, let us introduce the following 
loops on A D : 

S p (tt — s) — <£>(— s, ±</?o), sG[0,7r] for some (fo ^ 0(7r), (55) 

= *(-«, -2s), se[0,7r/2] . (56) 

S p is an oriented 5 1 -orbit, and S q is a loop which is everywhere transversal to 
the S 1 -action and oriented according to the flow of X 2 . 
The first quantisation condition: 



-!- / a + hn(8 p )/4E hZ + 0(h 2 ) 
27r Js„ 



is actually exact since Hi is a harmonic oscillator, and therefore reads: 

l-h = hn or h= 1 (57) 
4 ji + 3/2 

Because of the homogeneity property, we choose here to view this condition as 
a discrete quantisation of h (see remark |2.2| ). 

Assuming now that (j5^) holds, we can compute the semi-classical invariants 
associated to S q : 

The action integral is easily computed to be 

7T 

and the inva rian t e given by ( |l2] ) is equal to e/2. The sub-principal action I$ q 
of Definition [T^ is given by 

J<5 := / k = 3eln2. 



Finally, we can show by slightly shifting S q to the right (in the 9 direction) that 
its Maslov index is —2. Moreover 8 q turns around 70 in the direct sense, hence 
the regularised Maslov index is — 2 + (^ + n). 

We can now write down the second quantisation condition: Co = hol(5 q ), 
which reads: 



(1 + i{-l) n c-^) r Q + tej e £ (f 



(58) 



= e «(5fe+3eln2+f (_2+(f +n))+0(fc)). /gg) 

Using ( p7| ) and e = e/2, we obtain the equation in e and n: 
(1 + z(-l) n e - ^) r Q + e5(f- iln ("+ 3 / 2 )) = e <|e^2-fn+f+o(fc))_ ^ 
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Remark 3.2. The semi-classical invariants were computed explicitly; this 
is related to the fact that A Q admits a parameterisation as a "rational" vari- 
ety. Somewhat paradoxically, it would be more technical to compute the WKB 
invariants attached to regular tori, since no rational parameterisation of these 
tori exists. On the other hand, the regular invariants can be asymptotically 
recovered from the singular ones using Stirling's formula (J35|) . A 

Remark 3.3. The obtained formula (|60| ) yields easily (exactly as in pTfl ) 
the fact that the level spacings for the eigenvalues in a region of bounded e are 
of order 0(l/n\nn) — 0(h/ In h) - while they are of order 0(h) in a regular 
region. Moreover, the precise shape of the spacing function is readily derived 
and involves the log-derivative of the Gamma function. A 



3.2.5 Matrix form for H 2 

The goal here is to study the restriction of H 2 to the eigenspace <E n of Hi 
corresponding to the quantum number n (ie. to the eigenvalue h(n + 3/2)). 
In analogy with formulae ((5^), the operators l^j]) and ( |50|) can be written 

H\ = h (aAhMh) + 2a 2 (h)b 2 (h) - |) , (61) 



H 2 = V20 (a 2 {h)bi{h) 2 + ai (h) 2 b 2 (h)) 



(62) 



with ai(h) = {2h)- 1 / 2 (h^+x) and b x (h) = a x {h)* = {2h)- 1 ' 2 {~h^+x) (and 
similarly for a 2 (h) and b 2 (h) with the variable y). Using the unitary transform 
in L 2 (R 2 ) : f(x) — » ^/hf(Vhx), the operators Hi and H 2 are transformed 
into those given by the equations ( |6l| ) and ( |62| ) with dj{h) and bj(h) replaced 
by cij := Oj(l) and bj :— bj(l). Note that this shows that the homogeneity 
argument used for the classical analysis has an analogue in the quantum setting: 
if we know the spectrum for some value of h > 0, then the spectrum for any 
other value of h immediately follows. 

Now, using the Bargmann representation, we identify aj (respectively bj) 

k £ 

with the operator (resp. Zj), and let them act on the monomials which 

form a Hilbert basis of eigenvectors of Hi (corresponding to the eigenvalues 
Ei = h(k + 21 + 3/2)). Then it is easy to find the matrix of H 2 in this basis of 
£„ (n = k + 21) : 



H 



2 re, 



= V2hi 



( A nA 

A n , X X A nt2 

A n , 2 ^ _ 



(63) 



0/ 



with A n , t = v^(n -2i+l)(n-2£ + 2), I = 1, 2, . . . E[-]. 
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3.2.6 Numerical computations 

Since ft is of order 1/n, one sees that the coefficients of (|3|) are bounded as 
n — > oo. Moreover, since no coefficient A n j vanishes, the spectrum is simple. 
For these reasons, it is reasonable to expect a good accuracy of numerically 
computed eigenvalues. The resulting spectrum will be called the "quantum" 
spectrum. 

On the other hand, numerically solving equation ( |60| ) in the variable e - 
assuming that he remains in the bounds of the momentum map - yields the 
so-called "semi-classical" spectrum for Ei — 1 and h — (n + 3/2) -1 . If we wish 
now to fix ft and compute the rest of the joint spectrum, the same formulas ( pq ) 
and ( p^ ) can be used if one lets e = 2 J E and replaces h by h = h/E\. 

The results are displayed in the following figures. In Fig. [l^ we have su- 
perposed the quantum and the semi-classical joint spectra. The differences 
are hardly noticeable (they are theoretically of order ft 1 - which means ft 2 for 
the unsealed spectrum - and experimentally much better at the critical value 
H% = 0), even for a very large ft and very small E\% - both of these conditions 
are supposed to reach the limitations of our analysis. 



E 2 /h 




Figure 19: A comparison between semi-classical and quantum results. Here 
ft = 2/63 (so that Ei = 1 corresponds to n = 30) 

In the other figures we focus on one spectrum (here at E\ = 1) around the 
critical value Ei = - which is the most interesting feature. 
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;° (l//i)(E^(quantum) — (semi-classical)) 



b) 



(E. 



Ik) „(fc+l) 



)/2h 



Figure 20: The spectrum at Ei = 1, n — 30 (quantum and semi-classical are 
indistinguishable), a) The spectrum sorted in increasing order and displayed 
versus the eigenvalue number, b) The difference "quantum-semi-classical", c) 
The energy spacings. 



20 0° 3D 



Figure 21: Contrary to the numerical diagonalisation of the matrix (p3[), the 
semi-classical formula allows very small values of h — and the results are suppos- 
edly even more accurate. Here is displayed a window of size [— lOh, lOh] for the 
spectrum and the corresponding eigenvalue spacings, where we have let E\ = 1 
and n = 10 15 (h ~ 10~ 15 !! of course we haven't tried the matrix diagonalisation 
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3.3 Schrodinger Operators on S 2 
3.3.1 Setting of the problem 

We consider now the operator H = A + V where A is the canonical Laplacian 
on S 2 = {(x, y, z) £ M. 3 \x 2 + y 2 + z 2 — 1} whose spectral theory is given by the 
spherical harmonics and V : S 2 — > R is a smooth potential. We introduce the 
pseudo-differential operator H2 on S 2 which is obtained by averaging V using 
the 27r-periodic quantum unitary flow U(t) = exp(it^J A + 1/4) : 

Definition 3.1 

1 f 2 * 

H 2 = Tr u(-t)vu(t)dt . 

27r Jo 

The following results have been obtained by Weinstein and Guillemin (see , 
g§, 0, |§ and also g §): 

Theorem 3.2 • H2 commute with A. 

• H 2 is a PDO of order whose principal symbol is the Radon transform of 
V: 



H 2 (z) = — ^ V(^(z))dt , 

where ift is the geodesic flow with unit speed. The sub-principal symbol of 
H2 vanishes. 

• There exists an unitary FIO VL such that 

Or 1 HQ. = A + H 2 + R 
where R is a PDO of order —2. 

where the Hi 's are the orthogonal projections on the spaces Tii of spherical 
harmonics of degree I. 

Proofs can be found in @ and @. 

In such a way, we get a quantum integrable system H\ = ft 2 A, H2- The 
spectrum of H is related to the joint spectrum 

(h 2 l(l + 1), fii, m ), I = 0, • • • , 00, — I ^ m ^ / 

of (Hi,H 2 ) by Xi :Tn = 1(1 + 1) + fi^m + 0(l~ 2 ) and high energy asymptotics 
(I — > 00) corresponds in the usual way to semi-classical asymptotics h 2 . 1(1 + 1) = 
1, ft. — ► 0. We will study the system : 

Hi<p = h 2 l(l + l)ip — ip, H2^> = ehf , 

assuming that is a critical value of saddle type of H2- 

The Radon transform H2 of V is a function on the manifold Geod of oriented 
closed geodesies of S 2 . Geod is a global Poincare section for Hi and can be 
identified with S 2 C by associating to the circle t — > 7(f) = ucost + 

vsmt the unit vector u A v. Then reversing the orientation of 7 corresponds 
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to antipodal symmetry a on Geod — S 2 and H2 is even with respect to that 
symmetry. We can then interpret H2 as a function on the projective plane. 
This fact implies that, if H2 is a Morse-Bott function, it cannot have only 
local maxima and minima: it has always saddle points for which our analysis is 
needed. 

We will from now assume that we are in the simplest situation where H2 : 
Geod — > R has only 2 maxima, 2 minima and 2 nondegenerate saddle points. 
The singular manifold Ao is then the union of 2 tori which intersect along 2 
circles. The projection of Ao on Geod = S 2 , i.e. the reduction of Ao, is the 
union of 2 circles which are invariant by a and which intersect at 2 antipodal 
points. 




H 2 > 

Figure 22: G in Geod 
3.3.2 Semi-classical computations 

We will complete the computations in the simplest case where V itself is a 
generic harmonic polynomial of degree 2. Using 50(3) invariance, we have only 
to consider the 2-parameter family given by 

V a ,b,c{x, y, z) = 2(ax 2 + by 2 + cz 2 ) , 

with a + b + c = and a < c < b. 
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Because the Radon transform commutes with the SO(3) action, by Schur's 
lemma, H 2 itself is an harmonic polynomial of degree 2 on Geod which is given 

by 

H 2 = aX 2 + bY 2 + cZ 2 . 

The critical values of H 2 are a < c < b and H 2 — c = B 2 Y 2 — A 2 X 2 with 
A = \Jc — a and B = \Jb — c. It is clear that the singular value (1, c) of (H 1 ,H 2 ) 
is of hyperbolic type, so that we can apply the previous tools. 

Let us denote by G the projection of Ao on Geod by the map 7r which 
associates to a point m e U*S 2 the geodesic to which m belongs; G is the 
graph introduced in the general situation and is the union of the two circles 
C T = {AX = tBY, t = ±1}. We can compute the projections on G of the 
vector field X 2 and of the sub- principal form k on A (because k(A\) = 0). We 
will use spherical coordinates (6, ip), ^ 9 < 27r, ^ <p ^ 7r, on Geod: 

X = sin ip cos 9, Y = sin ip sin. 0, Z = cost/? . 

The symplectic form on Geod is 50(3) invariant and of total area Air. We will 
assume that: 

u> = sin ip d9 A dp . 
By direct computation we get on C T : 

d 

X 2 = 2tAB sin ip— 

dp 

and 

k = — — dip . 

2AB sin 99 

Let us denote by 7^, j = 1, • • • , 4, the four cycles of G oriented by X 2 , consisting 
of the union of one arc of C\ and one arc of C_i and bounding a topological 
disk in Geod. 71 and 72 are defined on figure ^ and 73 = (7(71), 74 = f(7 2 ). It 
is then easy to check using the explicit formula for the integral J dp/ sin^j that 



v;=v--,4,/ 73 . = ^io g -^. 



Moreover, we find easily, using formula (|l2|), that 



£2 = 2IB ■ ( 65 > 

We put a = atan-g and = atan^ = ^ — a. Action integrals are = 
A 3 = -Aj3 and A 2 = A 4 = 4a. 

Let Tj = hol(Tj-). We observe first, using the fact that h(l + i) = 1, the 
relations t\ — t 3 and t 2 = r 4 (modulo 0{h)). We get 

Tl = e i(-(4i+2)/3+/-Br+0(ft)) _ ^gg^ 

and 

T 2 = c i (( 4i +2)"+^+0(/i)) _ (gy) 

It follows that H := n = r 2 = r 3 = r 4 . It would be nice to prove that this 
relation holds mod 0(h°°). 
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Using the computations of |L0| p. 493 and putting T = T(e) with s = 
e 2 + 0(h) we get the following quantisation rule: 

de«( M -r(; J)r( J », V))=°- < 68 > 

Putting 



T = £ 



1 w 
cj 1 



we get: 



and the quantisation rule: 



H={1±lo) 2 £ 2 



j((4l+2) a +I+0(h)) _ J_tvI , • e(ir-2tlog(i+i)) ( 



n " r U + ie) W- AWg l ,+ a^ (1 ± ie~ £7V ) 2 , (69) 



which has to be solved in e where e enters in / and in e. 



Remark 3.4. Because of the ± sign in (|6C|), the spectrum can be separated in 
two spectra. For each of these, the spacings of eigenvalues are easily computed 
as in remark |3.3[ . Moreover, formula (^0|) shows that far from the critical value 
(ie. e — > ±00), the "+" and "— " eigenvalues associate in doublets, and that 
there is a universal transition happening when crossing the critical value (e = 0) , 
where a doublet "++" becomes a doublet " — h" . The details for these formulas 
are similar to 101 and left to the reader. A 



3.3.3 Matrix form for H 2 

Since V is a harmonic polynomial of degree 2, the Toeplitz operator H 2 is given 

by 

h 2 = n 2 yn 2 , 

where V here is just the multiplication by V and Tl 2 is the orthogonal projector 
on the space Tt 2 of spherical harmonics of degree 2. We shall first determine the 
explicit formula for H 2 with a generic V € Ti. 2 and then apply it to the specific 
form V = 2ax 2 + 2by 2 + 2cz 2 . 

The spaces Hi are seen as the spaces of the irreducible representation of 
so(3) on L 2 (S 2 ), acting via the differential operators 



" dz dy 

t V A tA. 

Lz — Z-q^ — XgJ, 

which are subject to the relation 

[L x ,L y ] — L z ( and cyclic permutations of (x, y, z) ). 

These operators commute with A and thus preserve Hi . As usual, (see eg. Q ) , 
we use the coordinates £ = x + iy and z, and let 

"± = iL x ± Ly, 
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so that tt + = -C|j + 2z|r and = fl* + = - 22^. A natural basis of Hi 
is then the following: 

<8, = (c',o + c',^C', . . • , ^C ; , • • • , n+'C 1 )- (70) 

We shall use the convenient equivalent representation given by the action of 
Slt(2) on the spaces V21 of homogeneous polynomials of degree 21 in C 2 = (£, 77). 
Using the following identification 



21 

' ()>l 



we get f2 + = £tt~i an( i a natural basis of V21 is the following: 



(^^ +?7 2 ',^,...^vV,-.-,W)- 

In the rest of the argument, this basis together with the basis (f70|) will be used 
to identify Hi and V21 ■ With this identification, V assumes the form: 

^=^(^ 4 +^)-3( a + 6)^ 2 . 

Up to a multiplicative constant (depending on I), there exists, for each Z, a 
unique equivariant morphism II : 7i 2 <£> 7~ti —*T~Li. Hence II is a multiple of T> 2 , 
where 

„ d d d d 

V = — <g> ® — . 

drj d£ <9£ drj 

(A subsequent composition by the multiplication f®g^ fg is always assumed.) 

On the other hand, using the fact that every element of the form fg, f G H2 
and g € TLi splits into 

fg = fi+2 + r 2 f l +r 4 f l _ 2 , 
where r is the radial distance and fj S H.j, one easily computes 

1 



/i = n(/® ff ) 



6 + 4Z 



- 2(2; - i) A « 3< ^ g ) 



Testing this formula and I? 2 with, for instance, f = x 2 — z 2 and g = (x + iy) , 
one gets 

n = V 2 

6(2Z- 1)(3 + 2Z) 



6(2Z - 1)(3 + 2Z) \d^®^f-~ 2 d^] ® d^] + lfrj 2 ®d£? / 



3 EXAMPLES 



58 



H- 



2[H, 



3(2/- l)(3 + 2/) 



(71) 



Now it is easy to let this operator act on VHi, and one finally obtains the 
matrix representation of H2 in the basis 25; 

N 

"S X "A 2l J 

with 

A m = 6(a+b)(3m(2l-m)- 1(21-1)) 
B m = 3(a-b)(m+l)(m + 2) 

Notice the symmetry of the matrix: A m — An_ m . 



3.3.4 Numerical computations 

As it is, the matrix (^3|) is very badly prepared for being numerically diag- 
onalised. Indeed, the spectrum exhibits near degeneracies - as is expected 
from tunneling effects - and usual algorithms will rapidly fail as / increases. 
Fortunately, there is an easy way to cope with this, for two commuting trans- 
formations can split the matrix: the projection onto the subspace spanned by 
the vectors from the basis ( |70| ) having even index, and the central symmetry 
(m — » 21 — m) of the matrix. We arrive at the following 4-blocks decomposition 
for the matrix —3(2/ — 1)(3 + 21)H2\h^ each block being a tridiagonal matrix: 



if I = 2k + 1 is odd: 

/ A a B a 

B2I-2 A-2 B2 

B21-4 \ 



''•A 2 k-2 B 2 fc_2 

-S21-+2 (A 2k +B 2k ) ) 



( A 

B21-: 



Bo 

A 2 B 2 

B 2 l~4 \ 





"■■P-2k-2 
B 2 k + 2 



B 2 k-2 

(A 2k -B 2 



I Ai Bi 

1 B 2l _ 3 A 3 B 3 



B21-5 \ 

' X B 2 k-3 

''■A 2 ) s _i 2B 2 fc-i 1 

B 2 fc+1 A 2 k + 1 / 



/ At Bi 

/ B 2l - 3 A3 



3^-83. 

B21-5 \ 



V J^2k-3 
B2k + 3 A2k-1 / 



if / = 2k is even: 



/ A B, 



B21-2 A 2 B 2 
B21 — A \ 

V B2k-& 
"'■■Alh — I 2B 2fc _ 2 

B2k A2k 



( A 

1 B21-2 A 2 B 2 



An B„ 

A 

B21-4 \ 




^ < E>2k-4 
B2k + 2 A2k — 2 / 



B2I- 



A 3 



\ ( 



''■A-2k-3 B 2 k-3 

B\k+i (A 2k -i+B2k-i) ) 



At Bi 

B 2 !-3 A3 B 3 .. 



(I 



'■A-2k — 3 B 2 fc_3 

6 2k + i (A 2 k-i-B 2 
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Figure 23: A comparison between semi-classical and quantum computations. 
Here is displayed the spectrum of Hi ^ in increasing order versus the eigenvalue 
number, for I = 40 and the potential V defined by a = — 1 and 6=1. The light 
crosses linked by line segments are the semi-classical computations while the 
quantum eigenvalues are the dark diamonds. We observe a very good accuracy 
(of mean order 0(h) = 0(1/1) predicted by the theory, but much better near 
the critical value). Notice also how the eigenvalue doublets reassociate when 
passing through the critical value. 



REFERENCES 



60 



References 

L.V. Ahlfors and S. Sario. Riemann Surfaces. Princeton University Press, 
1960. 

M. Audin. Courbes algebriques et systemes integrables : geodesiques des 
quadriques. Expositiones Math., 12:193-226, 1994. 

L.M. Bates and R.H. Cushman. Global aspects of classical integrable sys- 
tems. Birhaiiser, 1998. 

S. Bates and A. Weinstein. Lectures on the Geometry of Quantization, 
volume 8 of Berkeley Mathematics Lecture Notes. AMS, 1997. 

M.S. Child. Semiclassical Mechanics with Molecular Applications. Oxford 
University Press, 1991. 

Y. Colin de Verdiere. Sur le spectre des operateurs elliptiques a bicar- 
acteristiques toutes pcriodiqucs. Comment. Math. Helv., 54:508-522, 1979. 

Y. Colin dc Verdiere. Spectre conjoint d'operateurs pseudo-differentiels qui 
commutcnt II. Math. Z., 171:51-73, 1980. 

Y. Colin de Verdiere and B. Parisse. Equilibrc instable en regime semi- 
classique I : Concentration microlocale. Comm. Partial Differential Equa- 
tions, 19(9-10):1535-1563, 1994. 

Y. Colin de Verdiere and B. Parisse. Equilibrc instable en regime semi- 
classique II : Conditions dc Bohr-Sommerfeld. Ann. Inst. H. Poincare. 
Phys. Theor., 61(3):347-367, 1994. 

Y. Colin de Verdiere and B. Parisse. Singular Bohr-Sommerfeld rules. Com- 
mun. Math. Phys., 205:459-500, 1999. 

Y. Colin de Verdiere and J. Vey. Le lemme de Morse isochore. Topology, 
18:283-293, 1979. 

G. Darboux. Theorie generale des surfaces. Chelsea, 1972. 

J.J. Duistermaat. Oscillatory integrals, Lagrange immersions and unfold- 
ings of singularities. Comm. Pure Appl. Math., 27:207-281, 1974. 

A.T. Fomenko. Topological classification of integrable systems, volume 6 of 
Advances in soviet mathematics. AMS, 1991. 

V. Guillcmin. Some spectral results for the laplace operator with potential 
on the n-sphcre. Adv. in Math., 27:273-286, 1978. 

V. Guillcmin. Some spectral results on rank one symmetric spaces. Adv. 
in Math., 28:129-137, 1978. 

V. Guillcmin. Band asymptotics in two dimensions. Adv. in Math., 42:248- 
282, 1981. 

V. Guillcmin and D. Schaeffcr. On a certain class of Fuchsian partial 
differential equations. Duke Math. J., 44(1):157-199, 1977. 



REFERENCES 



61 



[19] F. Hirzcbruch. Topological Methods in Algebraic Geometry, volume 131 of 
Grundlehren der math. W. Springer, New York, 1966. 

[20] W. Klingcnbcrg. Riemannian Geometry, de Gruyter, 1982. 

[21] J. Moser. Geometry of quadrics and spectral theory. In The Chern Sym- 
posium, pages 147-188. Springer, 1980. 

[22] Z. Nguyen Tien. Singularities of integrable geodesic flows on multidimen- 
sional torus and sphere. J. Geom. Phys., 18:147-162, 1996. 

[23] Z. Nguyen Tien. Symplectic topology of integrable hamiltonian systems, I: 
Arnold-Liouville with singularities. Compositio Math., 101:179-215, 1996. 

[24] Z. Nguyen Tien, L.S. Polyakova, and E.N. Selianova. Topological classifica- 
tion of integrable geodesic flows on orientablc two-dimensional manifolds... 
Functional Anal. Appi, 27:186-196, 1993. 

[25] S. Vu Ngoc. Formes normales semi-classiques des systemes completement 
intcgrables au voisinage d'un point critique de Fapplication moment. 
Preprint Institut Fourier 377, http://www-fourier.ujf-grenoble.fr/ 
~svungoc/f n-tout .ps .gz, 1997. 

[26] S. Vu Ngoc. Sur le spectre des systemes completement integrables semi- 
classiques avec singularites. PhD thesis, Universite Grenoble 1, 1998. 

[27] S. Vu Ngoc. Bohr-Sommcrfcld conditions for integrable systems with criti- 
cal manifolds of focus-focus type. Comm. Pure Appl. Math., 53(2):143-217, 
2000. 

[28] A. Wcinstcin. Lectures on Symplectic Manifolds. Number 29 in Regional 
Conference Series in Mathematics. AMS, 1976. 

[29] A. Wcinstein. Asymptotics of eigenvalue clusters for the laplacian plus a 
potential. Duke Math. J., 44(4):883-892, 1977. 

[30] H. Weyl. The theory of groups and quantum mechanics. Dover, 1950. 
Translated from the (second) german edition. 



